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1.0 INTRODUCTION 


The purpose of this final report is to summarize the mathematical modeling 
and programming of the complete simulation program. Much of the material 
included is taken from earilier reports (Task 1 through Task U) but is 
repeated here in order to provide a complete description in a single doc- 
ument. 

Appendix A — Deck Setup and Input Variables - is also nearly identical to 
the separate document — Space Station Antenna Simulation Digital Computer 
Program Operating Requirements — and is included here both for completeness 
and to show the input parameter values used in the program run of Appendix 
C - Trade-off Study Simulation Results. 
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2.0 BASIC PROGRAM STRUCTURE 

The simulation program described in this report is an application of a gen- 
eral purpose simulation package (DEADAM) maintained by Dept. 62-11, LMSCI, 
Sunnyvale, Calif. This foundation program provides for reading input data 
cards, numerically integrating up to $0 first order differential equations 
(not necessarily linear or continuous), and monitoring up to U8 variables 
on printed output and on plots . 

For constructing a particular simulation program communication with the 
DEADAM package is accomplished through the subroutine SETUP in which the 
details of the differential equations and output variables are defined, 
usually with the aid of additional subroutines called from SETUP. 

The computations done by or called for by the SETUP subroutine for the "Com- 
puter Steered High Gain Antenna Simulation" are indicated in Figure A-l and 
are described in detail in the following sections of this report. 

For efficient computation of pointing requirements (gimbal angles), a simpli- 
fied simulation program is included which uses a different version of the 
SETUP subroutine ( SETUP/ PC)INT ) . The computations called for by this program 
are a subset of those illustrated in Figure A-l and are indicated in Figure 
A-2. 
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3*° SIMULATED SYSTEM MODELLING 


The system being simulated by the "Computer Steered High Gain Antenna 
Simulation" program consists of l) a high gain antenna, 2) an antenna 
gimbal control system, 3) an on board computer, and i*) the environment 
in which l), 2), and 3) are to operate. The interconnection of the 
mathematical models of these four subsystems is indicated in Figure 3.0-1 

3.1.1 Antenna Model 

This section describes the mathematical model used by the simulation pro- 
gram for the radiating system (antenna) and has been extracted from the 
Task 1 Report - Math Modeling of High Gain Antenna - dated 1 March 1972. 

The programming of this modeling is contained in the subroutines M0NELS 
and THERMO which are called by SETUP to computer gain along the actual 
line of sight, monopulse system output signals, and shift in the line of 
sight due to solar heating deformations j and in the subroutines GAIN, 
AXIS, and M0DEL which are called by the antenna initialization routine 
ANTIN only once (per case) to determine parameters used by M0NFLS and 
THERMO. 
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sensed S/S rate global angles & rates, S/S rates 



Figure 3.0-1 Interconnection of Subsystems 





3.1.2 AN OVERVIEW 

A simplified schematic diagram of the radiating system model is shown 
in Fig. A -2. 1. Each of the elements in this figure will be discussed separately 
in later sections. In this section we shall present an overview of the radiating 
system model and the interrelationship of the basic elements. 

The input portion of the program takes all of the input data needed for 
the calculation of gain except for the coordinate angles 0 and 0 and the coordinates 
of the sun with respect to the reflector. It has been assumed, in setting the input 
up this way, that the computer program will be used primarily in a "simulation" 
mode — that is, to simulate the operation of an actual or proposed system. This 
means that all the adjustable design parameters are introduced in the input block 
of the program and represent one particular system. All these input parameters 
will remain fixed throughout a run. The most rapidly varying coordinates then 
will be the LOS coordinate angles and the coordinates of the sun. Thus, with one 
input block of data we should be able to simulate the operation of the system 
throughout an orbit or an operational situation, but making parametric variations 
of certain design parameters (such as reflector diameter) for a fixed orientation 
of the LOS would require re-inputting the basic data. 

In addition to the basic design parameters, the input data base will 
contain a number of flags to select various options available within the program. 

One of the flags contained in the input data base is SMP which selects 
either the simple gain and beamwidth calculation or the more complex computation. 

If SMP=1, then the peak gain and beamwidth will be calculated from simple formulas. 
If SMP=0, then the peak gain of the antenna alorgits nominal boresight axis and the 
beamwidth are calculated using the current distribution method. The latter option 
will have several feed options and reflector options associated with it. 

i 

The output of either of these calculations is the peak gain of the antenna 
and the half-power beamwidth with the beam direction defined as being along the 
nominal boresight axis of the antenna. This output is then corrected for thermal 
effects. The angular position of the sun relative to the reflector is fed into a 
correction calculation block. Based on this sun position, the direction of the 
peak of the antenna beam will be shifted from its nominal position, using a table 
of values based on ATS F and G experience. The maximum gain value will also 
. be reduced in the same way. 
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INPUT 



Figure A-2,1 OVERVIEW OF RADIATINU SYSTEM MDDEL 












Following this correction, a pair of coordinate values 0 and 0 representing 
the LOS to the relay satellite will be entered. The computation of these coordinate 
angles will take into account any movement of the reflector boresight that results 
from steering — that is, the angles 0 and 0 represent the true position of the LOS 
with respect to the nominal boresight axis of the reflector. Using these coordinates, 
the gain in the direction 0 and 0 will be calculated using a simple approximation 
formula for the peak of the main beam. This gain is that associated with a single 
feed located at the focal point of the reflector. For either the non-monopulse 
option or the five -horn monopulse option, this would be the gain of the communications 
channel in the LOS direction. 

Another of the flags contained in the input data base is MON, which sets 
up the calculation of monopulse sum and difference signals where desired. For 
MONK), no monopulse signals are calculated and the gain value already determined 
is simply read out. Then the next values of 0 and 0 and sun angle are read and the 
calculation is continued. For MON=4, a four-horn monopulse system is assumed. 
Using design parameters contained in the input data base, sum and difference signal 
gain values are calculated. The sum channel gain value is also the communications 
channel gain value. After output, the program goes on to the next LOS position. 

For MON=5, the calculation is the same as for the 4 horn system, except that the 
communication channel gain value is the same as it would be for MONK). 

One very important point should be noted here. Even if the complicated 
form of gain calculation is selected, this computation is performed only once to 
determine the peak gain of the antenna and the half-power beamwidth. Once this 
computation is made, it is not repeated as new values of 0 and 0 are entered. The 
complex computation simply determines the maximum gain of the antenna and the 
relative gain at any angle is found using this peak gain value and an approximation 
formula. The gain at the actual position of the LOS could be determined directly 
from the complicated form of the computation, but the error in using an approxi- 
mation is comparable to (and probably less than) the computational accuracy of the 
complicated computation itself. The accuracy is certainly good enough for 
engineering purposes at any rate, and so it was not deemed prudent to use the 
complicated computation more than once per run as it requires a large amount 
of computer time compared with that required by the approximation. There arc 
some limitations to the foregoing statements, so we shall present the rationale 
here in detail. 
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3.1-3 ACCURACY OF THE GAUSSIAN PATTERN APPROXIMATION 


The current distribution method will be used (when desired) to calculate 
the peak gain on axis and the gain at four symmetrically located points off axis 
(in the vicinity of the half-power level). From these data an elliptical paraboloid 
will be fitted to the peak of the main beam. For the case of a symmetrical main 
beam aligned with the nominal boresight axis, the current distribution method 
will determine G m and , where 20 is the half-power beamwidth of the 
antenna pattern. Then the main beam within the half-power contour will be 
approximated by the formula 

G ( <*) = G m exp -p(a//9) 2 (3.1) 

where G ( a ) is the gain in a direction which makes an angle a with respect 
to the boresight axis. In the above p is a constant. 

The pattern shape given by Eq. (3. 1) is a Gaussian curve and this 
approximation is sometimes called the Gaussian approximation. When converted 
to decibels, the equation becomes that of a parabola. 

10 log 10 G (a) = 10 log 10 G m - 3.010 (« / f}) 2 (3.2) 

For this reason, the approximation is often called the parabolic approximation. 

For the more general case where the pattern beamwidths in two orthogonal 
planes are not equal, Eq. (3. 2) must be generalized to the form of an elliptical 
paraboloid. 

A study was made of the accuracy of Eq. (3. 2) compared to some 
typical theoretical pattern shapes. Specifically, the gain pattern functions 
chosen were 

g ( <*) = (? ni J - a -^ 

y u n 

where u =( ¥ D / X ) sin a 
the wavelength. These are the pattern shapes that result from circular apertures 
of uniform phase having an amplitude distribution of the form 

A (r) = (1 - r 2 )"- 1 V-' 


m 


(3.3) 


is the aperture diameter and A is 
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where r ($l) is the normalized radius of the aperture. 


When G{ a )/Gm as given by Eq. (3. 1) is matched to Eq. (3.3) at the 
peak (a = 0) and at the half power point ( a » )» the error curve shown in 

Figure A -3. 1 is the result. This is a plot of the equation 


m ' 

for n = 1 and n = 4. Although in fig. 3. 2 D/A was 200, there was no 
significant difference for 20 ^ D/A $ 500. From the figure it can be seen 
that the Gaussian approximation gives a gain not more than 0. 05 dB too low 
within the half-power contour, although it deteriorates rapidly when a > /9 . 


AG = -3.010 ( a / £) 2 - 10 iog 10 (I 


c 


(3.5) 


As we shall mention later, the computational accuracy of the current 
distribution method is about 0. 07 dB. While the Gaussian approximation would 
not properly represent specially shaped beams, such as those having flat tops, 
it should be an accurate representation of normal pencil beam patterns applicable 
to communications problems. The accuracy of this representation is better than 
0. 05 dB within the half-power contour and improves as we approach the peak of 
the beam. 

Beyond the half-power contour the Gaussian approximation begins to 
produce erroneous results. But for the computer-steered system, a misalignment 
of more than half a half-power beamwidth would cause serious doubts concerning 
the efficacy of the computer-steering method since the misalignment loss would 
exceed 3 dB. 

The Gaussian approximation will also be used in the monopulse case, 
not only for the data or communication channel, but also to form the monopulse 
error signals. Although initial misalignment during acquisition may be greater 
than fZ , giving some error in gain, the monopulse will track toward the beam 
peak to reduce the error to within acceptable limits when the acquisition phase 
ends and the self-tracking phase begins. 


The approximation will be used to save computer time since the current 
distribution integration process is quite costly compared to the simpler approxi- 
mation. The approximation is considered sufficiently accurate for practical 
engineering purposes. Beamwidth and the angle between the LOS and the beam 
axis will be read out to flag situations where the approximation begins to deteriorate. 
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Figure A-3.1 GAIN ERROR OF GAUSSIAN APPROXIMATION 





3-1-U INPUT 


The radiating system model receives data for calculation from four 
different sources: the built-in data files, external data files, the input subroutine, 
and other elements of the simulation program. 

The built-in data files are as follows: 

1. Standard Feed Gain Function. This is an experimental feed gain 
function data file which can be used in the program to calculate the reflector 
gain. This data file will represent a typical experimental feed designed to provide 
a tapered aperture illumination with approximately a -10 dB edge taper. 

2. Thermal Gain Degradation Data File. This is a data file representing 
the peak gain degradation of a typical unfurlable reflector antenna as a function of 
sun position. It is used, as explained in Section A -7, to provide correction for 
thermal distortion effects. 

3. Thermal Beamshift Data File. This is a data file representing the 
beamshift of a typical unfurlable reflector antenna as a function of sun position. 

It is used to provide for correction for thermal distortion effects. 

Any of the built-in data files may be updated as required to reflect 
a current design case. 

The one external data file is an experimental feed gain function data 
file which the user must provide. 

The only data that the model receives from other elements of the 
simulation program are the coordinates of the LOS and the coordinates of the 
sun in reflector reference coordinates. These data are the independent variables 
for each single run. When the model is run by itself without the balance of the 
simulation master program, these data must be supplied on cards. 

All of the remaining parameters necessary for the calculations 
performed by the radiating system model are supplied by the user on cards 
through the input subroutine. Some of the input parameters act as flags for 
selecting various options which may be exercised. The important input 
parameters are: 
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BMW. . . . This is the beamwidth factor. This is the product of the 

half-power beamwidth and the aperture size in wavelengths. 

See Eq. 5.3. 

CSN. . . .This parameter flags the feed option to be used, as follows: 

CSN ^ 0: Selects either the built-in or the external 

experimental feed data file. See Sections A -6. 3.4 
andA-6. 3.5. 

0 - CSN ± 999: Uses analytical function for feed gain function. 

Value of CSN is used in the analytical gain function. 

See Sections A -6. 3.2 and A -6. 3. 3. 

CSN 7 999: Uses function that provides uniform aperture 
illumination. See Section A -6. 3. 1. 

DBK. . . . This is the diameter of an area that blocks the reflector aperture. 

See Section A -6. 5. DBK is in inches. 

DIA. . . . This is the reflector diameter in inches. 

EFF. . . .This is the aperture efficiency in percent. 

FTV. . . . This is the feed to vertex distance in inches. If FTV=0, the 

feed is placed at the nominal focal point of the reflector. 

GFE. . . .This is normally the peak gain of the feed. GFE and CSN 

determine the exact feed option, as follows: 

CSN <-0; GFE <0: The built-in reference experimental feed 
pattern is used. See Section A -6. 3. 5. 

CSN < 0; GFE 2 0: The external experimental feed pattern 
is selected. See Section A-6. 3.4. 

0 ^ CSN ^ 999; GFE ■£. 0: Uses analytical function for feed 

gain function with peak gain value for feed computed 
internally. See Section A-6. 3. 2. 
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0 •£ CSN < 999; GFE > 0: Uses analytical function for feed 
gain function with GFE as the peak gain value for 
the feed. See Section A- 6 . 3. 3. 


GHZ. . ,. .This is the frequency in GigaHertz. 

HUB. . . . This is the radius of the hub of the reflector in inches. By 
hub is meant any portion of the reflector at the vertex which 
is either a hole in the surface or a surface area not meant to 
be part of the reflector (as, for example, a flat plate covering 
a hole). 

MON. . . . This is the monopulse flag, as follows: 

MON-O: does not calculate monopulse error signals. 

MON=4: Assumes a four-horn monopulse system. 

MON=5: Assumes a five-horn monopulse system. 

See Section 3*1* 


PAN. . . . Flags the choice surface. 

PAN •£(): Selects perfect paraboloid. See Section A- 6 . 4. 1. 

PAN > 0: Selects perfect Flex-Rib surface and sets the number 
of panels equal to PAN. See Section A- 6 . 4. 2. 

POL. . . . Selects polarization of the feed. 

POL = 1: dipole mode of linear polarization. 

POL = 2: slot mode of linear polarization. 

POL = 3 : right hand circular polarization. 

POL = h*. left hand circular polarization. 
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PPW. . . . This stands for points per wavelength and determines the 
rate at which the aperture is sampled to calculate the 
integral for gain. A built-in trap will require at least 
600 points to be used to evaluate the integral. See Section 
3.1. 


SMP. . . . This flags the simple or the complex model. 

SMP=0: The gain is determined by the current distribution 
method. See Section 3.1. 

SMP=1: The gain is determined by the simple model. See 
Section 3.1. 

TDT(6,3). . . . This is a data array which defines the method of 

extrapolating the gain degradation and beamshift data for 
thermal distortion to other reflector diameters and frequencies. 
See Section 3*1* 

XAP. . . . This is the cross over angle in the 0 plane for either monopulse 
case. 

XAT. . . .This is the cross over angle in the 6 plane for either monopulse 
case. 

XVP. . . .This is the cross over level in the 0 plane for either monopulse 
case. 

XVT. . . . This is the cross over level in the 0 plane for either monopulse 
case. 

(NOTE: At least one and no more than two of the parameters 
XAP, XAT, XVP, XVT will have a non-zero value whenever 
monopulse is called for. The monopulse may be specified by 
offset angles XAP and/or XAT or by cross over levels XVP 
and/or XVT (but not both). If only one non-zero value is 
given, the case is assumed to be symmetrical. Thus if 
XAP#0 and XAT=0 then the program sets XAT=XAP. Sec 
Section A -9. ) 
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3.1.5 PEAK GAIN AND BEAMW1DTH CALCULATION - SIMPLE CASE 


If SMP=1, the simple model is selected. This causes the peak gain to 
be calculated using the formula 



(5.1) 


where 1) = the overall antenna aperture efficiency, D = the aperture diameter 
and X = the wavelength. The wavelength in inches is calculated from the 
equation 


X = 11.802854/f 


(5. 2) 


where f = the frequency in GigaHertz. The quantities V , D , and f are 
represented by EFF, DLA, and GHZ, respectively, in the input subroutine. If 
EFF=0, the EFF is set equal to 0. 55 (considered to be a typical value). 

The half -power beamwidth @o is determined by 




— degrees 
D 


(5.:t) 


where ^ is the beamwidth factor represented by BMW in the input subroutine. 

If BMW=0, BMW is set equal to 70 (a typical value) in the simple case computation. 


Note that in the simple case the beam is assumed to be rotationally 
symmetric. 


There is no justification for the use of the simple case formulas given 
above on the basis of accuracy. These formulas have been accepted for years in 
the antenna design field as representing typical designs. Normally, they would 
be used because specific design information is lacking about the antenna to be 
used. On the other hand, the simple model will normally be quite satisfactory 
when parametric tradeoffs not related to the radiating system itself are being 
studied. 


By making the proper use of input data the simple model can be used to 
represent an actual antenna having known characteristics. Thus, if we have 
measured the gain and beamwidth of a specific antenna, we may use Eq. 5. 1 and 
Eq. 5. 3 to calculate 1) and £ (EFF and BMW) which may then be set ns input 
quantities. r 
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3.1.6 PEAK GAIN AND BEAMWIDTH CALCULATION - CURRENT DISTRIBUTION 

METHOD. 

3. 1.6.1 GENERAL EQUATIONS 

The gain of the antenna will be calculated by the "current distribution" 
method. The reflector is assumed to be illuminated by a point source feed (one 
producing spherical phase fronts at the reflector). The feed produces a surface 
current distribution on the reflector, and to find the total radiation at any point 
in the far-zone of the reflector, we sum (integrate) the vector contributions of 
elements of surface current. 

The current distribution method is a standard procedure in the antenna 
science. We shall omit the basic derivation which starts with the boundary 
conditions at the reflector surface, and Instead will use the radiation integral 
given by Silver* as our starting point. In the following the notation E will be 
used to indicate that E is a vector. 


The total field at a point (R, 0, 0) in the far-zone of the reflec tor is 


-3#CR 

E = e 


./Z 1 

+ 

V/M©M) e (9, 0) -j J T 

R 

2 n 

v- _ 


V X 


where tc =2 W / \ 

A - wavelength * 

?rad = *k e total radiated power 

fi and « are the constitutive parameters of the medium 

G^(0, 0) = the (scalar) gain function of the feed. 

e (0,0) is a unit vector everywhere defining the polarization of 
the feed pattern. 


*S. Silver (ed.). Microwave Antenna Theory and Design. McGraw-Hill Book 
Co. , New York, 1949, pl50. 
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The vector I<j> is the component of the complex vector I which is perpendicular to the 
direction of propagation. That is, 

L = I (i • I ) + 7 (I .~I ) (6.2) 

* 9 0 0 0 

The vector I is given by 

W |n x (I p x ;,] «" 3k<1 ' Ip ' lH) ^3 <«•»» 

The coordinate systems are shown in Fig. A -6. 1. The feed is located at the origin. 
Points on the reflector are located by spherical coordinates ( P ) or by 

rectangular coordinates (x, y, z). The variables of integration can be chosen 
from either of these two systems. The far-fieid observation point is located by 
the spherical coordinates (R, 0, 0 ). The corresponding rectangular coordinates 
(X, Y, Z) are also shown to aid in visualizing the transformations that follow. 

In Eq. 6. 3 I. is a unit vector in a direction toward a point on the 
reflector while is a unit vector in a direction toward the far field observation 

point. The quantities Gj and 6 have already been defined. The vector n is 
a unit vector normal to the reflector surface. 

In Eq. 6. 1 the product accounts for the back radiation of the 
feed in a direction toward the observation point. Normally in high gain antenna 
systems with well-designed, directive feeds, this contribution is negligible in 
comparison to the summed contributions of the reflector currents. Therefore, 
we shall neglect back radiation here. 


Neglecting back radiation, we find that the gain function of the reflector 
system is 


G (©, 0) = n 


— « 
J T 


( 6 . 4 ) 
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where rj is the radiation efficiency of the antenna. 
T 

*] r = P rad 


where P. is the total input power. The radiation efficiency accounts for 
only the dissipative losses in the antenna. 

The problem then reduces to the evaluation of the complex vector 
integral T 

3 . 1 . 6. 2 CALCULATION OF THE VECTOR I 

The vector I devaluated by numerical methods on the computer. The 
general expression for I is given by Eq. 6.3. The integration is performed over 
S . the surface of the reflector. Either spherical or rectangular coordinates 
may be used as variables of integration; we shall use the latter here. 

The quantity >/g in Eq. 6. 3 is numerically proportional to the electric: 
field strength of the feed Ss a function of direction. In this form there is an 
implication that only the main lobe of the feed radiation falls upon the reflector. 

The first factor in the integrand, jGf / p , is the feed pattern adjusted for the 

relative space attenuation to various points on the reflector. 

The bracketed expression in the integrand is the only vector part of the 
integral. This term accounts for the conditions of reflection at the reflector 
surface. The exponential term is the phase factor. For a perfect paraboloid 
the phase factor is reducible to exp ( _2 jkF )• a constant. 

The reflector influences the value of T through p , the distance from 
the feed phase center to a point on the surface, and through n , the vector 
normal. We shall discuss feed and reflector surface options in following 
sections and confine our attention here to the remaining parts of the integrand. 
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The key element to be considered is the vector V defined by 


V=nx(ipxe) 


( 6 . 6 ) 


into components parallel to the rectangular coordinate axes, the 
7 , and i„ are invariant in direction and magnitude as we 


If we resolve 
unit vectors i , 

move across tie reflector. Thus we may remove the aforementioned vectors from 
under the integral sign and (when we take real and imaginary parts) reduce the 
complex vector integral I to a set of six real integrals. Of course, we may resolve 
V into rectangular coordinate components regardless of whether we choose 
spherical or rectangular coordinates as the variables of integration. 


The vector e describes the polarization of the feed pattern as a function 
of direction. Since the reflector is in the far zone of the feed, the longitudinal 
aomponent of l e~ is zero. Even if this were not so, the longitudinal component of 
e would make no contribution to the integral since 


(ip 


)m 0 


(6.7) 


Thus, without any loss of generality we may take the polarization vector to be 


6 ~ i © B © + V 


( 6 . 8 ) 


with the constraint 


B. 


9 


+ 1 M : 


= 1 


(6.9) 
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We could have defined e in terms of its rectangular coordinate components. 
But then even the simplest polarization would require using all three rectangular 
components. We might be tempted to let e =1 (orT ) be a linear polarization, 
but this would mean that as the wave progresses towardme reflector, the electric 
field vector would remain parallel to the x (or y) axis and this is not typical of 
practical feeds. If we define e in terms of the spherical coordinate components 
of the ( P , & j i ) system, the longitudinal component may be discarded, but 
this definition has the effect of specifying radial and circumferential components 
of polarization on the reflector surface . 


Defining e in terms of the far-field spherical coordinates ( R , 0 , 4> ) 
accomplishes two things. First, the relationship between the feed polarization and 
the polarization of the resultant far-zone field of the reflector is obvious. Secondly, 
the 0 and 0 components of feed radiation are easy to visualize and are typical of 
practical antennas. 


When B = 1, Bv = 0, we have a linear polarization similar to that of a 

dipole antenna. We shall refer to this as the "dipole" mode of polarization. On 

the other hand when B^ = 1, B Q = 0, the linear polarization is similar to that of a 

slot antenna. Since B 0 and B^ can be complex and can each be functions of 0 and 0, 

it is possible to construct any arbitrary polarization vector. For example, circular 

polarizations may be formed by setting b = 1 /[2 and B / = + j /\f~2 

Q 0 ~ 


Although we originally set up e in terms of spherical coordinate components, 
we will eventually need to resolve e into rectangular coordinate components to 
perform the indicated vector operations. Also we need to find the 0 and 0 components 
of the far field radiation to compute . Thus we shall need unit vectors 

i R , i Q , and expressed in both reflector and far field coordinates. 

Letting 


^ Vy^ + 


( 6 . 10 ) 
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we have 


% = T P - T x * +i y If + *. f 


L = -I JL + I v - ^ + I XZ 

0 x y v z 


PA 


PA 


*0 


-i z + i y 

7 a 2 

A A 


and 


Ifl = ± x cos © - i sin 0 cos 0-1 sin 0 sin 0 


7 = -i sin 0 -i cos 0 cos 0 - i cob © sin 0 
i 9 X y z 




+i sin 0 -i cos 0 

y z 


Writing the vector n in the form 


n=i n + i n + i n 
xx 7 7 z z 
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( 6 . 11 ) 
( 6 . 12 ) 
(6. 13) 

<6. 14) 

(6.15) 

(6. 16) 

(6. 17) 



we have 


V = i 


+ i 


2 * < 8 

A 


B Q y -Bj _2L2 ) 
V P 


+ n 

_Z 

A 


z B n + xy B 


( z °9 


n 

_a 


( il) -Ji 

p A 


(y B e -xz I^_) 
P 


+ i. 


~ n x ( z B 0 + xy B qI ) - (A % ) 

p 


A P 

The surface element ds is 


d jL) 

p 


(6.18) 


dS = dx dy 


-n. 


( 6 . 19 ) 


Since the product Vds contains components of n in both the numerator and 
denominator, it is immaterial whether the normal is directed toward or away from 
the origin. 


For the phase factor we must evaluate the expression T where 


r = P (1 - ip . ljj) 


( 6 . 20 ) 


This turns out to be 

r = P - x cos 9 + y sin 9 cos 0 + Z sin 9 sin 0 


( 6 . 21 ) 
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Thus 


where V 


1 “:// 

is given by Eq. 


v a 

P n • 

z 

6, 18 and r 


-jk r 

dxdy 

by Eq. 6.21. 


( 6 . 22 ) 


We shall now consider feed and reflector options available in the model. 

3.1. 6 . 3 FEED OPTIONS 

The feed pattern is represented in Eq. 6 . 3 by the produc t\/G f ($ £ ) e ( £ ) 
The factor ^G^ represents the amplitude pattern (on a field strength basis) and 

the vector e defines the polarization. The phase factor for the feed is unity. 

The feed radiation at the reflector surface is represented by^G^ e e "^7 P 
which corresponds to a far-field representation of the feed radiation. The analysis 
assumes that the reflector surface is everywhere in the far zone of the feed. 

Various input parameters permit us to select from five different feed 
options. These range from a simple theoretical feed gain function to an exact 
description of the feed pattern which has been determined by measurement. The 
use of any of these options will be dictated by the amount of information that is 
available about the feed and by the purpose of making the particular computer 
run. 

3.1.6.3*! Uniform Illumination 

Uniform aperture illumination represents an important special case. 

This gives the upper bound of gain for in-phase aperture illumination and is the 
basis for defining 100 per cent aperture efficiency. There are two ways to 
calculate the maximum gain for uniform illumination. One is by the simple 
formula 


G 

m 



2U 


(6. 23) 
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which is simply Eq. 5. 1 with ri ±1 . This does not allow us to determine the losses 
due to aperture blockage or due to surface distortions. It assumes a perfect 
paraboloidal reflector of diameter D and an unblocked aperture. 

The second method is to let 

G f (^» O = c '. P( V*0 (6.24) 

in Eq. 6. 22, where C is a constant. When this equation is used, the current 
integration method may be used with any of the available reflector options. Also 
various aperture blockages may be used. Uniform illumination gives more weighting 
to the edge of the reflector and less to the reflector vertex than a tapered illumina- 
tion does. Since blockage is assumed to occur at the center of the aperture, using 
uniform illumination will tend to de-emphasize the effects of blockage. Since 
reflector options tend to have maximum deviation from a perfect paraboloid at the 
reflector edge, using uniform illumination will tend to emphasize the effects of 
surface distortions. 

Frankly, including the uniform illu min ation option in the current -distribution 
computation will be of limited value in the final program, but it is easy to include 
and may be of some value for diagnostic purposes. The use of uniform illumination 
with a perfect paraboloid and no blockage is an important reference case that can be 
used to check the computational accuracy of the computer program, for example. 

The result obtained can be compared with the value of obtained from Eq. 6. 23 
and the accuracy thereby determined. 

Thus the uniform illumination case will be included In the complex model. 
However, the gain function will not be that determined by Eq. 6.24 which would make 
the illumination uniform for the particular reflector model chosen. Instead the gain 
function will be that which provides uniform illumination for a perfect paraboloid 
with the feed at the focal point. For this case 


25 


LOCKHEED MISSILES & SPACE COMPANY 


or 


o f (0,0 


G f 



The Illumination function is then 



LS sec^ 
Dp 




which for a perfect paraboloid reduces to 



(6. 25) 


(6. 26) 


(6. 27) 


(6. 28) 


This option is called by setting CSN > 999 in the input subroutine. 
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3-1. 6. 3*2 Theoretical Feed Gain Functions 

An analytical function may be used to represent the feed gain function. A 
typical function is given by 


G f (0,*) = 2 (a + 1) cos m 0 (6.29) 

for 0 <0<'9O° and 

G f (0,4) =0 (6.30) 

for 9«lKl80° . 

Silver* uses this function in his treatment of reflector antenna properties. 

We have used this function extensively in deriving design curves for reflector 
systems. We have found that it gives good relative gain values when treating 
reflector antenna design problems parametrically. This would indicate that the 
pattern shapes are reasonable representations of typical feed patterns. Patterns 
for this gain function are shown in Fig. A -6. 2. 

The peak gain value obtained when using these functions is not correct. 

These functions produce gains which are higher than those obtained in practice. 

This is because a practical feed whose pattern shape matches the cos m 0 curve 
will have a lower gain than 2 (m + l) for two reasons. First the theoretical 
feed has no radiation in the rear hemisphere ( ip > 90° ) while a practical feed 
will have some backradiation, however small. Second, the theoretical feed has 
no forward hemisphere side lobe structure. The practical feed, however, may have 
sidelobes or distortions of the main lobe skirt which will lower efficiency even 
though the distortions or sidelobes do not fall upon the reflector. 

When the value of m is chosen to optimize the feed pattern for a particular 
reflector (a particular F/D ratio), the peak gain corresponds to an aperture efficiency 
of 83% or thereabouts for values of m up to 8. Typical practical values range 
from 55% to 65%. To be conservative the model will be set to adjust the gain of 

*Op. Cit. , p425 
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the feed function to the low end of the range. 

Eq. 6. 29 can be modified to be 

( 'P , 4 ) = L (m+l) cob 0 

3 

This option will be called by setting CSN in the range 0 C CSN ^ 999 in the input 
subroutine, in which case m will be set equal to CSN. The quantity GFE must 
also be set ^ 0. 

3. 1 . 6 . 3 . 3 Modified Theoretical Gain Function 

Another way of using the theoretical gain function and at the same time 
obtaining realistic gain values is to use the actual feed gain as determined by 
measurement. In this case 

) = Gp cos <// (6.32) 

0 V ^ 90° and 

( <M ) = o 

(6.33) 

& < 180°. The quantity Gf is the measured feed gain. 

This option may be called by setting 0 < CSN $ 999 and setting 

GFE > 0 . In the input subroutine GFE is the measured feed gain corresponding 

to Gp above. 




for 


for 90° < 


Since 55/83 is approximately 2/3, 

(6.31) 
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3. 1.6.3* U Experimental Feed Patterns 

We can improve the correlation between the predicted gain values which are 
computed by the program and experimental gain values that would be measured on 
the range by replacing the analytical function describing the feed gain function with 
the measured gain function of the feed. It is now relatively easy to determine the 
spherical coverage pattern of a feed in a form suitable for inclusion in the computer 
program. 

A first order approximation might be to fit a suitable surface to the two 
principal plane patterns of the feed. Although this is generally better than an 
analytical function, it is not quite as good as using the complete pattern contour 
of the feed sampled at close intervals. 

The experimental feed pattern data must be supplied to the computer in a 
format compatible with the program. The data file format to be described here is 
directly obtainable by measuring the feed pattern with Scientific Atlanta Radiation 
Distribution Equipment. If the feed pattern is measured another way, it will have 
to be converted to the standard format before it can be used with the computer 
program and no provisions for this conversion are included in the model. 

The coordinate system used in measuring the feed patterns is shown in 
Fig. A-6. 3. The feed axis is aligned with the ^ = 0 axis. The resulting coordinates 
(\Jj £ ) are the same coordinates defined in Fig. A-6. 1. We have developed equations 
for other orientations of the feed axis, biit so far we have never found occasion to 
use them. The difficulty in using the orientation shown in Fig. A-6. 3 is that the 
polarization rotates during the measurement of the pattern. Therefore, when the 
feed is nominally linearly polarized, the patterns are measured with a circularly- 
polarized transmitting antenna and a constant 3 dB polarization loss is assumed. 

When the feed is nominally circularly -polarized, the appropriate circular polariza- 
tion is used for the transmitting antenna. Then when the data file is used in the 
computer program, either dipole or slot polarization is used for the linear case 
or the appropriate circular polarization is used for the circularly-polarized 
case. 
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Figure A-6.3 FEED PATTERN COORDINATES 
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This procedure is admittedly an approximation with regard to the variation 
of polarization throughout the feed pattern. Equations have been developed for making 
a series of feed pattern measurements with different polarizations so that the vector 
e may be determined point by point throughout the pattern. This requires an added 
cost of measurement and an increased amount of computer time for each run and 
thus far such a complication has not been found necessary. 

With the set-up shown in Fig. A -6. 3 the Scientific Atlanta equipment samples 
the feed pattern every 2 degrees in 'p and in £ . Starting at ip - 0 and concluding 

at ip = 90 , the model tower is positioned in 2 degree increments. At each ip position 
a complete conical pattern is obtained by rotating in £ and sampling every 2 degrees. 
The data file therefore consists of 91 ^ scans each consisting of a start -of scan 
marker and 180 data points corresponding to the sampling in £ .At the end of 
the 91st scan, a start -of-scan marker followed by an end-cf-file marker occurs. 

The complete data file, therefore, consists of 16473 data points including markers. 

The data file as measured represents only the relative gain pattern of the 
feed. Each value in the data file must be corrected so that it represents the actual 
gain of the feed in that particular direction. To do this, it is necessary to make a 
separate experimental measurement of the peak gain of the feed. 

Having the complete feed pattern on a data file allows us to integrate the 
feed pattern to determine the feed directivity. This may then be compared with the 
measured gain value of the feed to determine the feed radiation efficiency. Also, 
it is possible to integrate the portion of the feed radiation intercepted by the 
reflector which may then be compared with the total radiation to determine the 
spillover efficiency. 

The data points representing feed radiation are quantized to the nearest 
quarter dB and are limited to within the range of approximately 0 to 40 dB. The 
quantization gives an error of at worst +0. 125 dB for any particular value of the 
feed pattern, but on the average the error is much less. 

Also, the feed pattern is sampled at 2 degree intervals in both angular 
coordinates and this does not agree with the sampling pattern used in the 
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integration process. Therefore, feed pattern values at the integration sampling 
points must be found by interpolation. Since the feed is generally a smooth, slowly 
varying function over the angle subtended by the reflector, the interpolation error 
is not significant. 

The accuracy of feed representation will be discussed further in a later 

section. 

The experimental feed pattern option may be called in the program by 
setting CSN < 0 and by setting GFE ^ 0 in the input subroutine. 


3. 1.6. 3 . $ Reference Experimental Feed Pattern 


This option is the same as the experimental feed pattern option, except 
that an experimental feed pattern data file which is considered to be typical of 
practical horn feeds will be built into the program. 

This option may be called by setting CSN < 0 and GFE < 0 in the input 
subroutine. 

3.1.6.U REFLECTOR OPTIONS 


The reflector is represented in Eq. 6. 3 by the scalar functiortfand by 
the vector o’ . The function p - f ( ip ^ £) describes the variation of the 
distance from the origin (feed) to a particular point on the reflector surface as 
a function of the variables of integration. If the surface is defined in the form 

* = *<*,y) then 


P = 



+ y 2 + (g (x, 



(6. 34) 


The vector n represents the variation of a unit vector normal to thSe 
reflector surface as a function of position on the reflector. The vector n can 
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be obtained by taking the gradient of any scalar equation equal to a constant which 
represents the reflector surface. Thus if 


H =P- f (0,1) = 0 

or if 

H = z - g (x, y) = 0 

(depending on how the surface is defined) 
then 


(6.35) 


(6. 36) 


n = V H (6.37) 

\VH\ 

Included in the radiating system model are two analytical representations 
of special reflector surfaces. One is the paraboloid and the other is the flex-rib 
reflector. (The latter is applicable not only to the flex-rib unfurlable reflector 
but also to various "umbrella” type reflectors). 

For any particular kind of antenna it is possible to develop a distortion 
theory embodied in an equation p = f ( p $ ) or z = g ( x,y) to 

represent the actual distorted surface in the presence of orbital thermal and/or 
dynamic effects. Then p ( or z ) and the vector n derived therefrom may 
be used directly in Eq. 6. 3 to calculate gain as it is degraded by the orbital environ- 
ment. This has been done for the flex-rib type of construction on the ATS F and G 
reflector project, for example. 

The difficulty is that the form of the p or z equation depends 
explicitly on the distortion theory formulated for each particular type of reflector. 
The only type of general treatment possible is to expand the reflector surface in 
a general power series or similar expansion, but this becomes cumbersome when 
one attempts to relate it to a specific design for which a distortion theory has 
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been or can be developed. 

Since the treatment of thermal effects has been limited in this study, we 
shall include only the perfect paraboloid and the perfect flex-rib surfaces. The 
effects of a thermally -distorted surface will be handled as we shall describe in 
a later section. The reader can be assured that if a particular reflector design 
is chosen and p or z equations can be derived representing the surface in 
the presence of the orbital environment, the basic model can be modified to 
include the effects of these distortions. It would also be necessary to modify 
the input subroutine to handle new input data. 

3 .I. 6 .I 4 .I The Perfect Paraboloid 

The equation for a perfect paraboloid of focal length F with vertex located 
at (x, y, z) coordinates of ( 0, 0, z D ) is 



(6.38) 


(6.39) 


Since z^ is a constant 


V H =Vz d = V(z + x 2 + v 2 ) 

4 F 


(6.40) 


and 


V H = i x + i v + i 
x 2 F y 2 F s 


(6.41) 
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and n can be found from Eq. 6. 37 with 



(6.42) 


This option can be called by setting PAN < 0 in the input routine. 

3.1.6.U*2 The Fsrfect Flex -Rib Reflector 

The Flex-Rib reflector may be called a sectoral approximation to a 
desired paraboloidal surface. The frontal view of this type of reflector is shown 
in Fig. A -6. 4. In each wedge-shaped sector the surface is approximated by a 
parabolic cylinder. 

The equation for a parabolic cylinder represented by the sector 
designated by s = 0 in Fig. A-6.4 is simply the equation of the x-z plane intercept, 
since the elements of this particular parabolic cylinder are parallel to the y axis. 
Thus for 6=0 

x 2 

2 “ Z D ~ — 

* P C 

where F c is the focal distance of the parabolic cylinder. 

If we rotate this panel through an angle £ = £ s 


(6.43) 


, we have the equation 


is. °°° + y 8lnC.-) 2 (6.44) 

panels then Eq. 6.44 becomes 

z = z D - ( x cos + y sin (6.16) 

4 F c 
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z - z 


D 


In particular if we have N 



8=5 


s = 3 
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and the s -th panel extends over the range 

(2a-l) — < | < (2s + 1) ~- 

We may define a parameter W as 

W = x cos 2JLS. jf y s i n SJL& 

■ N 


(6.46) 


This has special significance since it is the distance of a particular point on the 
reflector from the center of the aperture as measured, not radially, but along the 
centerline of its own particular panel. In terms of this parameter 


P = \/x 2 + y 2 + (8^- (6 . 47 ) 

V 4*c 

As in the previous section we may take z ^ as the constant defining 
the surface. Thus 


V H = V( z + — ) 


(6.48) 


and 


Vh - I JL cos 2jTja + I jl 2a + T 


2F, 


N 


y gp sin N 


(6.49) 


Then n can be found from Eq. 6. 37 with 


!VHj= _1\/w 2 + 4F 
2Fo V 


(6. 50) 
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The quantity is the focal distance of the parabolic cylinder. We 
must determine the relationship of Fq to F , the focal length of the "equivalent" 
paraboloid. To do this consider the intersection of a radial plane ( £ = a constant) 
with the reflector surface. Separating out the dependence on £ in Eq. 6. 45, 
we have 

r 2 2 ft 2trs. \ 

z = z n - 75 — cos - — B -) (6.51) 

D N 

where t - p sin 0 is the offset of a point from the z-axis. Comparing 
Eq. 6. 51 with Eq. 6.43, we conclude that the intersection of a radial plane containing 
the z axis and the reflector surface is a parabolic arc of focal length 

Fj (I ) = F C ( 6 . 52 ) 

cos 2 

N 

The focal length F^ of the reflector rib (which is the boundary between two 
adjacent panels) is obtained by setting ( = ( 2s ± l) 7r/N . Thus 



We may conclude that the paraboloid which is a best fit to this surface 
must intersect the surface twice in each panel and will therefore have a focal length 
: F such that Fr > F > F^ . As the number of panels increases 
without bound, the Flex -Rib surface more closely approaches a paraboloid. Since 
F r -> F c as N -=*• cd , then in the limit F R = F = F c 

A mean value of the function F^( £) can be found by integrating F^ over 
the aperture. There is, of course, no radial variation of . This mean value 
is 


F = F C (H-i- tsn 2 -£-) 
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or 



We have chosen to let F, the focal length of the "equivalent" paraboloid, be this mean 
value of the function ^ . 

When Zj) is set equal to F as defined above, the peak gain varies as a 
function of D/X as shown in Fig. A-6.5.When the reflector is a good approximation 
to a perfect paraboloid, setting Bq to any value other thah F does not improve the 
gain by a significant amount (something of the order of 10~ 3 dB). However, when the 
gain curve departs significantly from the perfect paraboloid gain curve, some of 
the lost gain can be recovered by moving the feed. This, however, is not a normal 
design region for this type of reflector since, even with refocussing, the gain 
degradation is severe. Thus, we can set z p = ^ and conclude that F is the 
focal length of the equivalent paraboloid for all normal reflector designs. 

The perfect Flex-Rib reflector case can be called by setting N > 0. 

The variable N is represented by PAN in the input routine. For this case p is 
found from Eq. 6.47 and the vector IT is found from Eq. 6.49 (normalized by 
Eq. 6. 50). The focal lengths are related by Eq. 6. 53 through Eq. 6. 55. 

3. 1.6.5 ASSUMPTIONS AND ACCURACY 

In this section we shall discuss the basic assumptions associated with 
using the current distribution method to calculate gain. We shall also mention 
some of the checks that have been made to determine the accuracy involved. 

There are three basic theoretical assumptions involved in the original 
formulation of the general equations. First of all, the method is applicable to 
determining the far field radiation of the reflector antenna only. Secondly, the 
reflector is assumed to be in the far field of the feed. Third, direct radiation 
of the feed to the point of observation in the far zone is neglected. 
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Figure A-6.5 PEAK GAIN OF FLEX RIB REFLECTOR 



The usual criterion for determining the far-zone of a directive antenna 
is that the distance R from the antenna exceeds 2D 2 / X . This will nearly always 
be true for orbital antennas. For example, at 100 nm altitude or at 100 nm distance 
between two space vehicles, it would require an antenna larger than 173 feet in 
diameter to fail to satisfy this criterion at X-band (10 GHz), clearly beyond the 
state of the art for the immediate future. The only conditions under which we 
might fail to satisfy the far-zone criterion are during the testing of orbital antennas 
on the ground on some ranges; during very close approaches of two space vehicles, 
as in rendezvouspr in some extremely large aperture systems intended for the 
transmission of power by microwave radiation. None of these situations are intended 
to be covered by the simulation of the antenna in this study. 


2 

When the distance from the antenna is exactly 2D / X , the gain loss is 
typically of the order of 0. 05 dB. Since normally the distance will be much in 
excess of this minimum distance, the gain error will be reduced to an insignificant 
value. 


The assumption that the reflector is everywhere in the far zone of the 
feed is not as easy to justify but normally it would be true. To illustrate, a typical 
horn type feed for a paraboloidal reflector would have an effective aperture of the 
order of 1 wavelength so as to illuminate the aperture properly. With F/D ratios 
of the order of 0. 5, this would mean that reflectors smaller than about 2 wavelengths 
would fail to meet the criterion. A 2X reflector fed by a lX feed would suffer 
severe blocking and provide only a few dB gain. In fact, it probably would be better 
to turn the feed around and eliminate the reflector, perhaps adding aperture to the 
feed horn. Thus, for normal feeds we will probably not run into a situation where 
the reflector is not in the far zone of the feed. 


The assumption that the reflector is in the far zone of the feed is made 
so we can apply geometrical optics when determining the amplitude of the feed 
radiation at the point where it strikes the reflector. Because of this assumption, 
there is no phase factor relating to the feed radiation in the basic formulation. 
The feed is assumed to be a point source with spherical phase fronts. There 
may be some cases where the feed radiation does not have equal phase across 
the angle subtended by the reflector and when this happens, the feed cannot be 
considered to have a single point phase center — thus the criterion would not be 
satisfied. 


U2 


LOCKHEED MISSILES & SPACE COMPANY 



Normally the criterion that the reflector be in the far zone of the feed will 
be satisfied to a degree sufficient for all practical purposes. The user should be 
on the lookout for any unusual feed designs which might have a significant phase 
error throughout the angle subtended by the reflector. Should this occur, he may 
correct the results for the appropriate amount of gain degradation attributable to 
such phase error or he may elect to use the simple model. 

The third assumption is that the back radiation of the feed can be neglected. 
If this is not so, the correction can be made as indicated in Eq. 6. 1. Normally, well 
designed feeds will have a front-to-back ratio of at least 15 or 20 dB and very often 
considerably more. For typical feeds , the back radiation would normally be some 
10 or more dB below isotropic and this would be insignificant in comparison with the 
peak gain of the reflector which normally would be 20 to 60 dB above isotropic. The 
back radiation usually would be further reduced in the normal design by obstructions 
behind the feed, such as feed support brackets, electronic boxes, or perhaps the 
space vehicle itself. 

Of more importance are the practical approximations which have been 
made in constructing the model from the basic formulation. For some of these 
approximations we have not been able to obtain an estimate of the error directly 
attributable thereto, but we know what the overall accuracy seems to be. 

We have made a basic assumption with regard to polarization. If the 
feed is linearly polarized, we have assumed that it is either a dipole or slot mode, 
as previously discussed. If the feed is circularly polarized, we assumed that 
the polarization is a quadrature combination of dipole and slot modes. No provision 
has been made for elliptically polarized feeds, since we normally try to achieve 
either linear or circular polarization. 

The choice of the polarization modes is only important when we are 
using the experimental feed patterns. We have made calculations using the 
theoretical feed gain functions, including uniform illumination, and have found no 
significant difference in the results for each of the four basic polarizations. 

When we are using experimental feed patterns, our assumptions regarding 
the polarization modes become of more importance. The difficulty is not whether 
the choice of dipole and slot modes is a correct one. We measure the gain function 
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of a linear antenna with a circularly polarized transmitting polarization and then 
during the integration process, we assume that the polarization is in one of the two 
standard linear modes. If instead the linear polarization has some other form, such 
as always parallel to the x-axis, we would still expect to get the same result, based 
on our experience with theoretical feed patterns referred to in the preceding paragraph. 
Similarly, if the feed polarization is circular throughout the pattern, it is immaterial 
what two linear polarizations are used to define circular polarization. 

If the feed has elliptical polarization, we shall be disregarding a component 
of its radiation either when we measure the feed pattern or when we use it in the 
program, or both. The error that occurs is dependent on relative amplitudes of 
the nominal polarization component and the omitted component. Since the reflector 
would have integrated the discarded component just as it does the retained component 
and would have then added the two in complex vector fashion, we can estimate bounds 
on the error from the relative amplitudes of the two components. 

Similarly, the four standard polarizations do not provide for the fact that 
some feed patterns have polarization characteristics which vary throughout the angle 
subtended by the reflector. This is most noticeable in frequency-independent type 
feeds over some portions of their operating bands. This variation tends to occur 
near the edge of the reflector, usually where the illumination taper tends to de -emphasize 
its effect. Furthermore, in a well-designed feed this variation is usually held within 
reasonable bounds. 

Frankly, it is difficult to determine what part of the error is attributable 
to the departure of the actual feed polarization from the assumed state. An investi- 
gation of this effect on accuracy would be a separate study all ip itself. We know 
that if it were a serious problem we could measure both amplitude and polarization 
(and even phase) of the feed pattern and put these corrections into the program. 

The result would be at least a four-fold increase in the test time required to measure 
feed characteristics, a significant increase in program storage requirements on the 
computer, and a significant increase in compute time for every case. The overall 
accuracy of the computation without this refinement is satisfactory for engineering 
purposes (as we shall see) so adding this complication does not seem justified. 
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The theoretical gain function option was primarily intended to be used as 
a means of evaluating parametric trade-offs. The relative accuracy for evaluating 
reflector designs is quite good and has been confirmed by experiment. In this test 
the gains of a high-precision 4 foot paraboloid and a 12 panel, 4 foot Flex -Rib 
reflector were compared at 4, 8, and 12 GHz. A direct substitution method was 
employed so that the same aperture illumination, focal distance, and frequency 
were used for both cases. The primary feed patterns were measured and fitted 
to the analytical feed functions. Then the relative gain values predicted by theory 
and experiment were compared. Agreement was within about 0. 05 to 0. 1 dB, 
even though the feed patterns were not exact fits to the analytical functions. This 
test confirmed the relative accuracy of the current integration method of gain 
calculation. 

As mentioned before, use of the analytic function without modifying the 
peak gain value gives peak gain values that are unrealistic. The modification given 
in Eq. 6. 31 should be more representative of typical antennas and should be within 
+ 1 dB for typical cases, though it may be in error by more than that amount for 
any particular case. When the analytic function is modified by using the actual peak 
gain of the feed, the error is estimated to be within + 1 for the particular case 
considered. In these estimates, we have not included any allowance for transmission 
lines losses and have referenced the antenna gain to the feed terminals. The estimates 
in this paragraph relate to absolute accuracy — that is, the accuracy of predicting 
the true gain of the antenna. 

The basic radiation integral is approximated by numerical integration 
methods in the model. The sampling interval must be fine enough to sample accurately 
the illumination taper across the aperture, the reflector contour, and the electro- 
magnetic wavelength. Electromagnetic theory would tell us that there should be one 
sampling point per half-wavelength across the reflector, but if this is used as a 
criterion, the reflector surface may not be adequately sampled when the reflector 
is small in terms of wavelengths. Thus a wavelength criterion is used, but a certain 
minimum number of points must be taken for the double summation to represent the 
integral. 
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For large reflectors, fewer points can be taken if the surface is smooth. 

The paraboloid, for example, requires fewer points than a many-panelled Flex -Rib 
because the paraboloid is a smooth surface of revolution. In many cases when doing 
parametric studies we have found good relative accuracy with sampling rates as low 
as one point for every five wavelengths, but, of course, the absolute accuracy 
deteriorates somewhat. 

The criterion we normally use is that one point should be taken per 
wavelength, except that in no case should fewer than about 600 points be taken. 

The accuracy for one point per wavelength sampling is almost as good as that 
obtained for one point per half-wavelength sampling, but the integration time is 
only one quarter as long. For exploratory parametric studies, fewer sampling 
points can be taken to reduce compute time while searching for the optimum values. 
The proper minimum sampling rate for a given problem can be determined by 
increasing the sampling rate and noting where the gain value begins to stabilize. 

There will be some oscillation of the gain value about the true gain value as the 
sampling rate is increased further due to basic inaccuracies in the calculation 
(computer arithmetic accuracy). 

To check the computational accuracy, we have used the uniform illumina- 
tion feed option and the perfect paraboloid reflector option and compared the resultant 
gain value with that obtained by using Eq. 6. 23. For one point per wavelength 
sampling the agreement is about 0. 07 dB. This is considered to be the computational 
accuracy of the program. (The machine used for this test is the SDS 940 time-sharing 
computer . The case considered was the ATS F and G reflector at 8.25 GHz). 

When using the experimental feed option, two additional errors occur. 

First the feed gain values are quantized to the nearest quarter-dB. Second, since 
the sampling of the feed pattern in the data file is not coincident with the sampling 
of the aperture in the integration process, feed pattern data must be interpolated 
at the integration sampling points. Studies of the best method of interpolation 
showed that linear interpolation was adequate since the feed pattern generally is 
quite smooth for typical designs. Moreover, taking higher order differences 
eventually and/or occasionally results in poorer accuracy due to the quantization 
of the feed pattern data. The reason for this can be seen by comparing the following 
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difference tables for unquantized and quantized data. 

Actual Values Quantized 

4.5 4.5 


4.6 

.1 




4.5, 

.0 




4.7 

.1 

.0 



4.75 

.25 

.25 



4.8 

.1 

.0 

.0 


4.75 

.0 

-.25 

-.50 


4.9 

.1 

.0 

.0 

.0 

5.0 

.25 

.25 

.50 

1.00 

5.0 

.1 

.0 

.0 

.0 .0 

5.0 

.0 

-.25 

-.50 

-1.00 


Due to the rounding off, higher order differences eventually will increase causing 
poorer accuracy of interpolation. If we used only first and second order differences 
we would occasionally get poorer accuracy (near where the values change). 

To check the overall computational accuracy of the program, we constructed 
an artificial "experimental feed" data file simulating the feed gain function that provides 
uniform illumination of the aperture (Eq. 6. 25). The feed data file was quantized to 
the nearest quarter-dB, as an experimental data file would be. In the integration 
process linear interpolation was used. At 8. 25 GHz the gain program calculated a 
gain of 57. 873 dB for a 30-foot reflector using this data file with sampling at one 
point per wavelength. The gain equation for uniform illumination (Eq. 6. 23) gave 
a gain of 57.958 dB. The difference, 0.085 dB, represents the overall computational 
accuracy, including machine accuracy, the representation of the double integral 
by the double summation, linear interpolation of the feed data, and the effects of 
quantizing the feed data. In this check a 500 panel Flex-Rib reflector was used to 
simulate a paraboloid, so the minor effect of that approximation is also included. 

Thus the computational, accuracy using the experimental feed option is almost the 
same as that obtained using the theoretical feed options. 

Aperture blockage by the feed assembly or other structures is accounted 
for by not integrating the surface currents over the projection of the blockage on 
the surface of the reflector. That is, the limits of integration are adjusted so as 
to omit portions of the reflector radiation not contributing to the on-axis gain. 

This approach is accurate on axis and near the refelector axis. For observation 
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points at wide angles to the reflector axis, it is necessary to account for the fact 
that the blockage would not obscure the vertex area of the reflector but some other 
displaced part. However, since the model is normally used to predict gain values 
near the axis of high gain systems , the refinement applicable to wide angles has 
not been included. 

Blockage is considered to be a circular blockage centered about the 
reflector axis. If the blockage is not circular, the best procedure is to equate the 
area represented by the blockage to a circular blockage area for use in the program. 
The blockage may not be centrally located, but normally it would be near the center 
of the aperture. If the blockage is off center, then the model tends to give conserva- 
tive results. 

The apparent gain of the reflector antenna will be modified by the impedance 
mismatch of the feed when mounted in the reflector. This mismatch may be thought 
of as having two components, one the inherent mismatch of the feed in free space 
and the other the reaction of the reflector on the feed. The theoretical gain of the 
reflector/feed can be corrected if the mismatch of the feed in the reflector is known. 
Often this is not available. Therefore, the model calculates the magnitude and 
phase of the interaction of the feed and the reflector. This can be used as a bound 
or tolerance on the gain values obtained from the model. Normally, in a well 
designed system this will not be a problem, since the feed should be reasonably 
well matched in free space and since, the normal aperture distributions for maximum 
gain do not provide a high degree of reflector feed interaction. 

We shall cite an example to illustrate the reaction effect and indicate, 
at the same time, something of the overall absolute accuracy of the model. LMSC 
experienced a perturbation in a gain versus frequency curve of about 1. 5 dB for two 
measurements 50 MHz apart at L-Band. Careful tests indicated this perturbation 
was not due to measurement error. The model was used to predict the gain at the 
two frequencies in question from experimentally determined feed patterns. The 
computed values were within 0. 1 dB of the curve obtained by extrapolating smooth 
portions of the gain curve with no correction for mismatch. This represented an 
agreement of about 0. 7 dB between theoretical and experimental values. When the 
reflection coefficient was calculated, it was found to be of the order of -11 dB with 
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a large phase change over the 50 MHz frequency interval. This reaction in combination 
with a measured free space mismatch of the feed of about 2:1 was sufficient to account 
for the observed gain perturbation. These were range instrumentation antennas with 
-22 dB edge illumination for low sidelobes. The low edge illumination resulted in 
relatively high concentration of energy near the center of the reflector causing more 
interaction with the feed than would normally be encountered in a typical general 
purpose antenna. The typical case where the illumination tapers to -10 or -11 dB 
to maximize gain results in a reflection coefficient attributable to the reflector of 
the order of -18 to -20 dB. 

One effect which is not taken into account is the interaction of the feed 
support struts with the feed. This is the effect the structure has in modifying the 
amplitude or phase of the feed pattern (as distinguished from the blockage of the 
aperture caused by the feed and other support structures). Normally, the antenna 
designer will design the structure for minimum effect on the feed patterns, but 
structural considerations nearly always limit the designers ability to achieve his 
desires. This is probably the major factor limiting the accuracy of the model. 

With well-designed and simple support structures we may neglect this consideration. 
But the user should be on guard when complicated support structures are used (as 
when large masses are supported with the feed). Where this effect is not negligible, 
the feed pattern can be measured in the presence of the support structure to obtain 
a better feed model. 

We have previously cited tests of the relative accuracy of the model 
(the tests with the 4-foot reflectors) and of the absolute accuracy of the model 
(tests on the range instrumentation antennas). We have also checked the absolute 
accuracy of the model using a feed designed for the ATS F and G verification 
tests and a four-foot reflector. The agreement between the gain value predicted 
by the model and that obtained by measurement was about 0. 16 dB. The conditions 
of test were such that we could expect very good experimental accuracy, so this 
test seems to confirm the accuracy of the model. 
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3-1.7 


CORRECTION FOR THERMAL EFFECTS 

The scope of the effort in considering influence factors was limited in 
this study to be compatible with available funding. As a consequence, the effects 
of thermal distortions on the gain of the antenna along a particular LOS will be 
handled as a correction based on extrapolation of LMSC's experience with the 
ATS F and G reflector and with other systems. 

The effect of a distorted reflector can be accurately handled by including 
in the current distribution integration model an option for a distorted reflector. 

The overall program sequence can be modified to calculate the actual gain of the 
reflector along the line of sight, rather than calculating peak gain and using the 
Gaussian approximation as described in Section A-3. All that is necessary to 
accomplish this is an analytical function which describes the distorted surface 
and which can be differentiated to determine the normal surface vector. 

The problem of including the precise calculation of distorted surface 
effects is the complete lack of generality in any formulation that can be made. 

The description of the perfect paraboloid may apply to many types of reflector 
construction, such as honeycomb sandwiches or spun dishes. Similarly , the 
description of the perfect Flex -Rib surface also applies to the umbre 11a type of 
deployable reflector. But to include distortion effects requires developing a 
distortion theory for each possible method of realizing each type of surface. For 
example, the distortion equation of a Flex-Rib reflector would be different from 
that of an umbrella reflector which might have rigid ribs that fold up against the 
reflector axis. Similarly, there is a large amount of work involved in establishing 
the design parameters that optimize each type of reflector for the chosen orbit. 

Thus in a limited scope effort, only a correction-type treatment of thermal effects 
is possible. 

The ATS F and G reflector designed for 8. 25 GHz operation (30 foot 
diameter aperture) will be the primary initial basis for making a thermal 
correction. Beamshift and gain degradation data for this case will be extrapolated 
to other sizes and frequencies as required. 

The peak gain and the on -axis gain of the ATS reflector at 8. 25 GHz is 
shown as a function of orbit hour in Fig. A4-7. 1. The beamshift is shown in 
Fig. A-7.2. Both these curves are for the configuration utilizing Chromel-R 
mesh. The orbit hour is directly relatable to sun angle in the xz plane of the 
reflector. (In far-zone coordinates this would be the ZY plane or the 0 = + if piane. 
See Fig. A -6. 1). In the simple or the complicated gain calculation, the peak gain is 
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Figure A-7.1 GAIN DEGRADATION - ATS F&G REFLECTOR 




8.25 GHz Chromel-R mesh 



Figure A-7.2 BEAMSHIFT - ATS F&G REFLECTOR 





calculated and the beam position is defined as 0 = 90 , 0 = 90°. The thermal 
correction will correct the peak gain value and provide for a new reference 
position for the beam. A coordinate transformation (a rotation about the reflector 
axis) will be used to handle the case where the sun is not in the aforementioned 
plane. 

We shall denote the gain degradation for the reference case by A G 0 for 


diameter D Q and frequency f Q . 
use the form 


Then for any diameter D and frequency f, we will 


10 log 10 G m = 10 log 10 G m ~A G 


dB 


(7.1) 


where G mo is the peak gain for the undistorted case and where 


A G = A G n 




D 2 

P 

o 


][ 


+ \ 5 


+ b„ 




(7. 2) 


In Eq. 7. 2 the parameters and b n are input constants which fit the design case 
best. The form of Eq. 7. 2 allows some flexibility of adjustment by making changes 
in the input constants. Generally, some of the input constants may be zero. The 
form used above does require some justification. 

Since most of the systematic distortion errors are monotonic, increasing 
functions of radius, we would expect the surface error to increase as the diameter 
is increased. As the surface error increases, so does the gain degradation. Thus 
we can expect the magnitude of the gain degradation to increase as D/D 0 increases. 

Eq. 7. 2 gives us the possibility of making a second-degree polynomial fit to the 
scaling effect. Similarly, frequency changes act in the same way because wavelength , 
which is inversely proportional to frequency, is the unit of measure in reflector 
antenna systems. (That is, all distances may be measured in Wavelengths instead 
of inches). 

To use the correction option, we must simply make some judgement 
about how the distortions of the case in hand would be extrapolated from the 
reference case and then assign appropriate constants. If the gain degradation 
data for the case at hand becomes available in some way, the reference file may 
be updated and then with an appropriate choice of constants no correction will be made. 
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The gain degradation will depend on other factors besides D and f. 
Amplitude distribution, focal length-to-diamter ratio, and various other parameters 
may all have a bearing on the actual gain degradation. But our intent here is to 
supply the best correction that can be made for thermal distortions consistent with 
the scope of the study. The manner in which this correction is made allows the 
program to be updated whenever distortion effects are analyzed in more detail. 


Exactly the same philosophy is applied to the beamshift problem. When 
the sun is in the xz plane (Fig. A -6. 1), the angular corrections to the position of 
the peak of the beam become 


6 = f + aA 

IT V 

and 


(7.3) 



e = 

n 




(7.4) 


where 



d 0 



(7.6) 


Using the same functional form for beamshift does not necessarily imply 
that the beamshift will be a function of diameter scaling ratio or of frequency scaling 
ratio. For example, a distortion mode corresponding to a simple tilt of the refloptor 
if scaled to a larger diameter, would still be a tilt of the same angle and thus would 
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be independent of D/D 0 . Also, we would not expect the displacement of the beam 
to be dependent on frequency when this displacement is measured in degrees or 
radians. That is because we tend to think of the displacement as the result of an 
"effective" rotation of the reflector. The importance of this rotation arises when 
we measure the rotation in beamwidths and not degrees. Thus while we have 
included provision for dependence on diameter and frequency, there may be many 
cases in which the beamshift would be Independent of either variable or both. 

Tim reference file for beamshift data can be updated just as we mentioned 
above for the gain degradation reference file. 
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For purposes of providing the input signal to the monopulse model (subroutine 
MgfMFLS) a distorted antenna reference frame (D-frame) is defined as the frame 
obtained from the undistorted antenna reference frame (A-frame) by applying 
the (small) rotation dtp and d0 about the appropriate axes (see definitions of 
dip & d0 above). In this reference frame the peak gain of the antenna lies 
along the — Z D axis and the unit line of sight vector in this reference frame 
(LHATD) is the input to M0NFLS. 


The D-frame is obtained from the A-frame by an effective single rotation 


d t =V(<*P) 2 + 


2 - - A dco - 

(d0) about the 6^ axis (see Figure A.O.l) where ^ € 


d6 M 

+ ^0* The Eu ^ er parameters for the transformation from the A-frame to 

the D-frame are defined such that 


L q 3 


p- “I - dir _ df 

= hl A siji t * %- cos 


or since dt is small and e^ is a linear combination of £ and €g. 


V 

q 2 

Si™ 

II 

cos 0 
sin 0 

♦# 

- sin 0 
cos Cl 

1 

O' | 


0 


0 

- _ 


i (dtp cos 0 — d6 sin Q ) 
q_ = i (ckp sin Cl + d0 cos Cl ) 


q^ = 0 


q u 


= 1 


whs re Cl is the angle between the axis and the axis of the d < p rotation 
(see Figure A.O.l). The transformation matrix is obtained ;frora these 

Euler parameters using the subroutine EFT0DC . 
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3 . 1.8 


GAIN CALCULATION ALONG THE LOS 


The first step after the input subroutine is to calculate the antenna gain 
at five points. The pattern of these five points is shown in Fig. A-8. 1. When the 
current distribution method is used to calculate the gain* the gain is calculated on 
axis and at four angles oft axis as shown in the figure. When the simple method is 
used, we have the same result since we know the peak gain on axis and the half- 
power beamwidth. Thus we have the following table of values in either case: 


Gain 

0 

0 

G, 

0 

0 

1 

n 

n 

g 2 

»n 

in *8 

G 3 

®n 

+ 8 

°4 

®n “8 

«n 

G 5 

®n + 8 



In the above table, ® n ,0 n are the coordinates of the reference axis (in this case the 
beam peak). When the method has been used to calculate gain, G 2 , G 3 , G 4 , and 
G5 have the value Gj-3. dB and 8 is the half-power beamwidth divided by 2. 

The gain values and angular coordinates may have been corrected for thermal 
effects. 


These five gain values and their corresponding coordinates are used to 
fit an elliptic paraboloid to the peak of the beam. This paraboloid is then used to 
determine the gain along the LOS for varying coordinates of the LOS, thus saving 
much computer run time. The accuracy of this approximation has been discussed 
in Section 3. 

In the vicinity of 0 = 7T/2, 0 = 7T/2, Eq. 3. 2 may be generalized to 

the form! 


10 l Ogl0 G a - 10 log^Q G (0, P) = T (0-<y 2 + P(P -0 O ) 2 (8.1) 

where G m is the peak gain value and where 0 O , 0 Q are the coordinates of the 
beam peak. It will happen that G m - Gp 0 Q = 6 n , and 0 O = 0 n because the 
reference axis 0 n , 0 n happens to be the beam peak in these calculations. 
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Figure A-S.1 LOCATION OF GAIN CALCULATION POINTS 
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Iii general, however, it is not necessary that this be so, and when the current 
distribution is applied to distorted reflectors, the peak gain and the beam position 
can be found by using Eq. 8. 1 and five gain values as indicated in the above table. 

For convenience we shall use the notation G' x to indicate 10 log 1Q G v 
in the following. 

In general, we find that 


P= it 2 (2G * ‘ G ' 2 “°3 j 


( 8 . 2 ) 


T = 


26 


<*4 * G i 


-g 5 ; 


(8.3) 


and that 


= J»n + 


°3 ~ G 2 


4P4 


(8.4) 


G 5 " G i 

0=0+ -a. 

o n 4T , 


(8.5) 


When the beam peak lies along the reference axis then G *2 = G' s and G' 4 = G' g 
so that 


10 


lo «io 0 <•» W - <£ - - (G^-y ^e_-_9 o ,J 


0 I - P - T 

60 


( 8 . 6 ) 
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S B* 


which for the calculations made in the simple model would be 

10 log^ 0 (9, pi = 10 log^, Gj - iai2 [(W 0 ,! * < 8 - *</] < 8 - 7 > 

AM* 

where 2 is the half-power beamwidth. Equation 8. 6 is used to represent the 
antenna pattern in all subsequent calculations. The parameter 8 is chosen 

the current distribution calculation to be half the estimated half-power beamwidth 
obtain good correlation over this region* 
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3.1.9 


MONOPULSE MODEL 


The monopulse system used on the high gain antenna could be either a 
four-horn or a five-horn configuration. In the five-horn configuration a single 
horn located at the focus of the reflector provides the communications channel 
signal; we shall denote this horn and its associated beam by the letter T. The 
error signals are formed from four offset feeds a, b, c, and d clustered about 
the focus. In the four-horn configuration the on-axis horn is omitted and the 
remaining horns are not displaced quite as far from the reflector axis. In this 
case the communications signal is the sum signal of the monopulse. 

The choice between the two configurations is a tradeoff between 
communications channel gain and tracking capability. The sum signal of a four-horn 
monopulse, being the sum of four offset beams, has a lower peak gain that that 
provided by a single on-axis horn. Thus the five-horn configuration provides 
higher communications channel peak gain. On the other hand, because of the 
presence of the fifth horn, the remaining four horns must be displaced further 
from the axis of the reflector than they would be in the four-horn configuration. 

As a result, the four-horn approach can have higher beam crossover levels — 
more nearly at the optimum level — and provides better tracking. 

The monopulse model will calculate four gain values for each position 
of the LOS. One of these, G D* is the gain available along the line of sight in 
the communications or service channel. Another, Gg, is the gain available in 
the monopulse sum channel. The remaining gain values, Gq and Gg, are the 
monopulse error channel gains for the constant 8 and constant 0 "planes" 
respectively. The model allows the user to select a four -horn, or five-horn 
monopulse or no monopulse at all. The flag for this is MON in the input routine. 

When M0N=0, the model will not compute monopulse sum or difference 
signals. Then Gjj is set equal to G(8,0) as discussed in the previous section. 

The remaining gain values Gg, G 0 , and Gg are set to zero. 

When MON=5, G^ is also set equal to G(8,0). The model then 
calculates Gg, Gq, and Gg as we shall discuss in the following paragraphs. 

When MON =4, the model computes Gg, Gq, and Gg and then sets G D = G g . 

Each of the five beams can be represented by an equation of the 
form of Eq. 8. 6, except that 9 0 , 0 O must be adjusted to include the shift of 
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each beam off the axis represented by 0 n , 0 n . The coordinates of the axes of the 
five beams in the model are 


Beams 


0 


o 


t 

a 

b 

c 

d 


0- - U 

n 

9 n - U 
0 n + U 

0 + U 

n 


0 - V 

*n + v 

K - v 

*n * v 


For the small offset angle used in a monopulse we would not expect the 
peak gain G j or the beamwidth to change. 


(NOTE: In the previous presentation of the monopulse model in the 
progress reports a rotationally symmetric beam was assumed similar to Eq. 3.3. 
It was then necessary to find the angle between the LOS and each of the 5 beam 
axes. The approximation represented by Eq. 8. 6 gives the gain in terms of 
0, 0 of the LOS directly. In the vicinity of 0 = 90°, 0 = 90° we can treat 0 and 
0 much the same as we would treat rectangular coordinates. Thus with Eq. 8. 6 
as the model the geometrical problem is simplified) 


In the direction 0, 0 we obtain G a , G^, G and G^. Since sum and 
difference signals must be obtained by adding and subtracting on a voltage basis, 
we have for the sum channel gain 

°s= A a^r®r + Vv + ^r + v°r ) 2 <9 - i> 

where A g is a factor to account for the modification of the peak gain by the 
monopulse action. 

Similarly 


G 0 = A 

0 s 


+ Vs - v°c~ “V^y 


(9. 2) 
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v \ (V°7 + V^ -V°4 ) a <! 

v 

Quadrant can be determined by the phases of the difference signals compared to 
the sum (on a field strength basis). 

If we let 

- PV 2 - TU 2 
lOw 

where P and T are given by Eq. 8. 2 and 8. 3 and where 


W = log l6 e “ 0. 434294482 (9.5) 

and if we let 


0 = 0 (9.6) 


] 


( 9 . 4 ) 



0=0 + 40 

n 

then 


(9.7) 


G = A G exp 
s 8 p r 



( A 


£][ 


4 cosh 


PVA0 

lOw 


cosh 



2 

( 9 . 8 ) 
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A G exp 
s p 


- P ( ? ff* co8h Pvi i elnh < 9 ' 9 > 

L 1Q W JL J ,? 


TUA0 

lOw 
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G 0 = A s G p ex P " P (A 0) 2 + T ( A 0) 2 1 (*4 sinh PVA0 cosh TUA9 ~j (9. 

*■ lOw -* L lOw lOw J 


10 ) 


and 


G t = G 1 . exp 



(9.11) 


The only actual gain value needed is Gjj. The tracking system will 
utilize the error signals normalized with respect to the monopulse sum signal. 
Thus the output of the monopulse model will be Gjj, G Q /Gg, and G^/Gg (plus 
the quadrant information). The required error signal ratios are 



(9.12) 


(9. 13) 


For the 4-horn case the sum channel of the monopulse is the service 
channel. For this case it is necessary to evaluate Ag so that Gjj( 9,0) = Gg(9,0) 
may be found. For the no mdnopulse-or the 5 horn-monopulse we do not need to 
evaluate Ag. 
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The quantity Ag may be found by performing the integration 



G g (0, 0) sin 9 d 9 d 0 = 


4 rr 


(9.15) 


Unfortunately, the model for Gg(6,0) is not valid for all 0, 0. But 
the model is best near 0 n , 0 n where its contribution to the integral is most 
important. A reasonably good value of Ag can be obtained by limiting the range 
of integration to a small region around 9 n , 0 n which we shall denote by r . 
Specifically, 

4 v 


16 A G 
s p 



(Q, 0) 
16 K 

s p 


sin 9 a 9 d 0 


(9.16) 


We already know G. , the maximum value of G T (0, 0) as given by 
Eq. 9. 11. We can find the value G- that the above process would yield. Thus, 

y 


G = 

y 


4 IT 


fir*- 


(o, 0) 


sin 0 d 0 d 0 
4 it 


(9.17) 


JJ exp P (A j) + T (A9)^ j sin 0 d 9 d 0 

The ratio of Eq. 9. 16 to Eq. 9. 17 is the gain loss due to the monopulse action 
as compared to that that would be provided by a single horn on axis. 


Thus, the maximum value of Gg (0 n , 0 ) is 


G (9,0) = 
s ' n n 


( 16 A b % > 


n 


(9.18) 
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If the model were perfect, then Gy would be equal to and G g (0 n , 0 n ) would be 
the value obtained from Eq. 9. 16. By using the ratio 16 Ag G p /G to modify Gj_, 
we obtain a better model since some of the inaccuracies in Eq. 9.16 and Eq. 9. 17 
tend to cancel. 

Finding the peak value of Gg is only a problem in the four-horn case. 
When MON=4, the integrations indicated by Eq. 9.16 and Eq. 9.17 will be 
performed for the first set of values 0, 0. The peak value thus obtained for 
Gg (0 n , 0 n ) from Eq. 9.18 will be retained for all subsequent LOS positions 
for the same case. 

3*2.1 Antenna Control System Model 

This section describes the mathematical model used by the simulation pro- 
gram for the antenna control system and has been extracted from the Task 1 
Report - Math Modeling of High Gain Antenna - dated 1 March 1972. 

The programming of this modeling is contained in subroutines CNTRL7 and 
SERV0. 
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3. 2. 2.0 


MODES OF OPERATION 


3 . 2.2.1 REFERENCE FRAMES 

In order to discuss the several modes of operation to be simulated, 
three of the orthogonal, righthand reference frames that will be used are described. 
The antenna reference frame (A) has, as two of its axes, the negative boresight axis 
and the inner most gimbal axis. The IMU reference frame is fixed with respect 
to the space station. The gimbal angles are Euler angles relating the angular 
position of the antenna reference frame to the IMU reference frame, assuming 
rigidity of all members. The outer gimbal angle (about the axis common to 
the boom and gimbal ring) is a and the inner gimbal angle is /? . The IMU 
senses the angular position and angular velocity of the IMU reference frame with 
respect to an inertial reference frame (I). 

>.2.2.2 COMPUTER COMMAND MODE 

In the computer command mode the antenna angular position command 
is computed by the on-board digital processor as a function of the ephemeris of 
the relay satellite at which the antenna is to point and the position and attitude 
of the space station with respect to the inertial reference frame as determined 
by the inertial measurement unit (IMU) and the on-board space station ephemeris. 

A spiral search scan may be superimposed for acquisition of the relay satellite 
RF signal. For slewing to a relay satellite, limits will be placed on commanded 
angular rates commensurate with mechanical limitations. The antenna angular 
position command is converted from components in the IMU reference frame 
to gimbal (Euler) angles prior to being compared with the actual antenna position, 
as shown in Fig. B-l. 


The rate stabilization shown as a block in Fig. B-l and in detail 
in Fig. B-2 provides the rate feedback path. In order to reduce the disturbance 
input effect of space station attitude changes on the antenna angular position 
with respect to the inertial reference frame, the space station angular rate is 
added to the antenna angular rate in the IMU reference frame. 
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Figure B-l MANUAL AND COMPUTER COMMAND MODES 
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Figure B-2 RATE STABILIZATION 






3. 2. 2. 3 MANUAL COMMAND MODE 

It is assumed that the manual mode is used to bring the antenna 
to a fixed orientation relative to the inertial reference frame and that the 
manual command signals originate within the space station by someone 
making appropriate settings of input devices and then initiating motion. 

Only step angular position commands relative to the inertial reference 
frame will be simulated. The commands will be converted to gimbal angles 
and then compared with the sensed gimbal angles to form the position error 
signal as shown in Fig. B-l. 

3. 2. 2.1* RF AUTOTRACK MODE 

In the RF autotrack mode, using a sufficiently strong RF signal 
being received from the relay satellite, the angular position difference between 
the line-of-sight to the relay satellite and the antenna boresight axis is sensed 
by the monopulse system. The servo is automatically switched from the search 
scan computer command mode to the RF autotrack mode when the received 
signal strength exceeds a preset threshold level. Since the receiving antenna 
beam shape changes as a function of carrier frequency, a corresponding error 
signal gain adjustment must be made as indicated in Fig. B-3. Next, the error 
signals are converted from the antenna reference frame (the frame in which 
they are sensed) to the gimbal (Euler) angles. 

Whenia the RF autotrack mode, the gimbal angle sensor outputs 
may be used to correct the on-board program which computes the gimbal angle 
commands required to point the antenna at the relay satellite, i.e. , to update 
the computer command mode program. 
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Figure B-3 R-F AUTOTRACK MODE 












3. 2. 3.0 


CONTROL LAWS 


3. 2.3.1 OBSERVATION OF OUTPUT STATES 

The control laws (more than one will be modeled and simulated) 
will use the angular position and rate of the antenna about the gimbal axes when 
in the computer command and manual command modes. When in the RF auto- 
track mode, the angular position error signal is developed within the monopulse 
system but the gimbal rates may still be used directly. Fig. B-4 indicates 
several alternate sources of information for the control law. When in the 
computer command mode, gimbal axis angular rate information may be sensed 
by a rate gyro or a tachometer, or estimated from the gimbal angular position 
sensor output. 

3. 2. 3. 2 STRUCTURE 

Control law alternatives are shown in Fig. 3 $ and include integral- 
error compensation, derivative-error compensation, and derivative -output 
compensation. The control law outputs, which are the compensated error 
signals for each gimbal axis, are the inputs to the motor drive electronics. 

As is indicated in the figure, the math model includes both linear and non- 
linear control law alternatives. The control law selected is implemented 
within the three "gain and compensation" blocks shown in FigB. B-l, B-2, 
and B-3. 

3. 2. 3. 3 VISIBILITY 

Possible obscuration of the relay satellite by the earth or by 
portions of the space station is monitored using visibility discretes. If the 
antenna servo system is operating in the monopulse autotrack mode and the 
line-of-sight to the relay satellite is obscured, it will be possible to 
automatically switch to the computer command mode. 
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Figure B-5 CONTROL LAW ALTERNATIVES - SINGLE AXIS 








3.2.U 


DIGITAL PROCESSOR 


The sampled-data effects of digital processing will not be simulated 
because the Adams -Bashforth integration package which will be used, has a 
variable integration step size. Simulation of a fixed sampling rate controller 
would introduce a discontinuity at each sampling instant which would require a 
re-initialization of the Adams -Bashforth integration package, thus substantially 
increasing the computer run time. This also precludes simulation of random 
effects using a random number generator. The entire servo simulation will 
be deterministic. All digital processing will be modeled as continuous signal 
processing. 
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3.3*1 On Board Computer Model 

This section describes the mathematical model used by the simulation pro- 
gram for the on board computer and has been extracted from portions of the 
Task 2 Report — Mathematical Model for the Pointing Requirements of the 
Space Station High Gain Antenna — dated 28 Oct. 1972, and the Task h Report 
Math Modeling of Interface Computer — dated 30 March 1972. 

The programming of this modeling is conained in the subroutines 0BKNMT, 
0BEIH3, 0BEFKL, fifeEIH2, INTSS1, HJTSS2, 0BSTCL, AHTCQM, and CLBRTN which 
are called by SETUP to compute the antenna gimbal angle commands predicted 
by the an board computer. 
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3*3*2 Computed Space Station Kinematics 

The on board kinematics computations are simulated by the subroutine 0BKNMT . 
Tito function of these computations is to determine the space station attitude 
relative to the orbital reference frame from the initial attitude and the 
time history of the space station and orbital reference frame inertial rates. 

Within the subroutine the relative space station attitude is specified by the 
euler parameters (q^, q^) of the rotation which would align the orbital 

reference frame (0) with the space station reference frame (c). When this 

A 

rotation is by the angle <p about the axis €, the euler parameters are defined 
as 

■ € x sin (<p/2) 

■ € ■ sin ((p/2) 

" sin ((p/2) 

= cos ((p/2) 


where 

€ 

x 

€ 

y 

6 

l z 


A A 

- [€] = [€] 
o c 


q l 

*2 

q 3 

q U 


78 

LOCKHEED MISSILES & SPACE COMPANY 



The tine derivatives of these euler parameters may be expressed in terms of 
the components of the space station and orbital reference frame angular 
velocities (oi and U) respectively) as 


• 

*1 

• 



*2 

• 

■ 5<°c +5 o> 

q 2 

q 3 

• 


q 3 

X 


%_ 


where 
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and the direction cosine matrix for the orbit frame to platform frame trans- 
formation may be expressed in terms of the euler parameters as 
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C Vo ] 


q l *4 ^3 2(q l «2 + q 3 \ ) 2(q l q 3 - q 2 \ } 


2(q x q 2 - q 3 q u ) -q 3 - q 2 -q^ 

J(q l q 3 * q 2 2 < q 2 q 3 “ q i 


2(q 2 q 3 + q x q^) 

2 2 2 . 2v 

" q l "*2 ^3 *V 


These formulations are derived from the classical works of Euler, Rodriques, 
Hamilton and others and are documented in a Lockheed internal report by 
E. F. Whiten. 1 


The above differential equations for the four euler parameters are numerically 

integrated and used as indicated above to obtain the direction cosine matrix 

[T > ] as a function of time, 
c/o 


1 


A Four - Vector Method of Describing the Rotational Motion of Rigid Bodies , 
E. F. Whiten, Feb. 19Sii. 
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3.3.3 C ctnputed Ephemerides 

The relay satellite ephemerides are computed by subroutines 0BEPH1 & 0BEHI2, 
and the space station ephemeris is computed by subroutine 0BEHi3« The com- 
putations are identical in all three of these routines. 

The on-board ephemeris subroutine 0BEPH1 computes the position vector R and 
the inertial velocity vector 7 by means of the approximate formulae that 
would be used by an on-board computer. While there are several possibilities 
for on board ephemeris formulas this first version of 0BEPH1 ( 0BEPH1/LAGRAN ) 
uses 8 point lagrangian interpolation for each of the components of R and V. 


For each of the 8 time points, values are supplied at initialization for the 
time (t^), the X, Y, Z components of R (Y^, Y^, Y^ respectively), and the 
X, Y, Z components of V ( Y^ j , Y^, Y^g respectively). 


The components of R and 7 are then computed as functions of time (t) according 
to 


[R] E - 

1 

- ^ - 

z u 

z 5 


L z 3 J 


h- 


8 


A 8 * ‘ t ( i 


l - S Q Y . Q = i 

i J ji* 77 t -t 

31 j=l 1 J 




For the space station the orbital reference rate (tD,) and the earth frame to 
orbit frame transformation matrix (CTq^J) are computed from 

5. = -ri- (R x 7) 

R 
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This expression for U is an approximation which neglects the very small com- 
ponent along the orbital z-axis due to gravitational acceleration along the 
orbital y-axis (an earth oblateness effect). The expression for [T Q y^] in 

terms of R and V is however exact, with U) being used as an intermediate 
. r 


function of R and V. 
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The object (relay satellite or ground station) located at R is considered 
■ ■ a 3 

to be visible from the station located at R if the line of sight between 
f s 

the station and the object does not pass through the sphere of radius 
REQ + 1^0000 ft centered at the center of the earth. This model neglects 
the oblate shape of the earth (a variation of about 70000 ft) and assumes 
no effect upon the line of sight above an altitude of lJJOOOO ft. 

The compulations are performed as follows: the angle (0! o ) between R g and R^ 

is compared with the angle ( 0 ^) defined by 
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and if 0^ > then the object must be visible and no further computations 

are required. If a <0; then the object will not be visible if it is on the 

o m 

opposite side of the horizon from the station, i.e., if the projection of 

R on R is less than P 
os o 


R . 
o 


P = 
o 


R 


R 


< P 


(REQ + 1^0000) sin a = ^■ REQ +1 ^9 00Q ) I 

m 


R 


£ . £ < (REQ + 150000)' 
o s 


The condition for line of sight intersection with the earth (approximation) 
is then 


a < a and R r R < (REQ + 150000)' 
o m os 
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3*3*5 Space Station Occultation of L0S 




Discussion 

The objective of this function is to determine whether the line of sight 

(specified by the unit vector S) orginating from the antenna gimbal point 

(specified by the position vector R) passes through the space station 

2 2 2 2 2 2 

ellipsoid envelope (specified by the relation x /a + y A + z A = l). 

The method used is to analytically determine the minimum value of the 
quantity 

_ A 2 / 2 , 2*2 2, 2 

Q = x /a + y A + z A 


along the line 


R + S d 

x x 


y • = R + S d 

■ y y 

z = R + S d 
z z 

(where d is the forward distance along the line of sight from the antenna 
gimbal point), and then to observe whether or not the minimum is less than 1 
and occurs at a positive value of d. Both of these conditions being met 
indicates an intersection. 

Setting the derivative of Q with respect to d equal to zero indicates that 
Q is minimum when d has the value 


S R /a 2 + S R A 2 + S R A 2 
x x y Y z z 

(S x /a) 2 + (Ta? + (S a A)* 


and the minimum value of Q is 


2 2 2 

( R + S d \ /R + S d \ /R + S d \ 

-j 2 - j 


es 



In terms of these quantities the consition for line of sight intersection 
with the space station envelope is ' 

d > 0 and Q <1 
o o 
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3.3*6 Gimbal Angles and Rates (ANTC0M) . 

Discussion 

The equations which form the basis for ANTC0M are given below. The gimbal 
angles are obtained from the line of sight (L0S) vector coordinatized in the 
earth frame (E) assuming a roll, pitch (R/P) rotation from the platform frame 
(C) to the antenna frame (H). The object is to find the gimbal angles (a and 
(3) which cause the L0S to align with the negative z axis of the antenna frame. 


L = RR — RA = L0S vector coordinatized in E frame 
E E E 

e. - < 1 , l>£ Magnitude of the L0S vector 
EE 


L = T_ L = L0S vector coordinatized in C frame 
C C/E E 


Tr/c = = Transformation from C frame to H frame 


H/C 

Ca] 


m 2 - 


*1 

0 

0 1 

0 

cot 

sa 

0 

» 

sot 

ca 

C/3 

0 

-S/3 

0 

1 

0 

S/3 

0 

C/3 


^ a rotation by a about the x axis / Sot ^ sin 

(. A 


•GOi = cos 


intt\ 
os a' 


= a rotation by /3 about the y axis 


L 

H 


Defining L 
C 


= L0S vector in H frame after the gimbal rotations 


*1 


u 2 


A 
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’ o" 


"c/3 S/3 Sa -S/3 Ca“ 


V 

0 

= 

0 Ca Sa 


i 2 

4 


S/S -c/3 sa C/3 ca 


u 


Solving the above equations for a and /S results in 



(Ca) i 2 = -(Sa) i 3 



Where the signs are chosen to align the negative 
z u axis with the L0S vector for /3 = 0. 

ii 

Similarly 

(C/3) + (Sa S/3) l 2 - (S/3 ca) = 0 

and = [(Sa) /? 2 - (ca) tan /? 

Noting that (Sa)/ 2 - (ca) = (Sa) (Ca) !£ 2 - (ca) 2 3 

ca 

(sa) 2 ^ 3 - (ca) 2 / 3 . / 3 

ca ca 
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Hence 




where the signs were chosen (for Oc = 0) to align the negative z u axis with 

n 

the L0S vector. 

The determination of signs in this equation can also be arrived at from con 
sideration of the third equation: 

- t = . S0 - C0 SGr + JL C0 C<* 


The expressions for the gimbal angle rates (a, 0) can be derived from the 
following starting point. 

Assume that the following variables are known 


C— I 

U) = angular rate of the platform frame (C) with respect to an 
inertial frame (I) 

I 

RR = time derivative of the relay orbit radius vector with respect 
to an inertial frame 


I 

RA = time derivative of the platform orbit radius vector with 
respect to an inertial frame 



transformation from earth frame (E) to plaform frame (C) 


0c> t h/c 


as defined previously 
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Since 


then 

and 


L = 
I 

L = 
C 

L = 
H 

L = 


C 

= L 


(RR-RA) 

( 1 1 \ 

\RR — RA/ 

-(a, 0 ' 1 x l) + L 
l + x , l) 

- (c/-° x l) 


It is desirable to coordinatize the above equation in the antenna frame (H) 
H ■ ♦ 

since both L and L have simple representations in that frame. Assuming the 

antenna to be tracking the L0S perfectly 

l = range 


L = 
H 


0 

0 

-/ 


H 

L = 
H 


0 

0 

l 


t = 


range rate 


Hence 


H 

L 

H 


L - x L\ 

H H H # 


where 


C 

L 

H 


3j 


(for convenience) 
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H-C 
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Cp 0 — S>3 

0 1 0 


» • « 
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<• a. 


S/S 0 cp 
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»- «<■ 
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D -S/3 a p 
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Cp a 
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Combining the above 

’ 0“ 
0 

-/ 

■> «< 

• • 

Solving for a ana /3 

P = 

a = 



+~2 
jPcjS 



range rate ^ 



3.3*7 Calibration Model (CLBRTN) /. 

The calibration model is a "best fit" of thB error between the open loop 
generated gimbal angles and those gimble angle profiles which exist when the 
antenna is tracking via the mcnopulse RF link. The term "best fit" will be 
used in the sense of least square. The procedure for obtaining the best fit 
is described in the sequel. 


When the antenna is tracking via the RF monopulse link the actual gimbal 
angles will be used as a reference. If the computer generated gimble angle 
commands were exactly correct there would be no difference between the gen- 
erated command and the actual gimbal angles of the autotrack mode. However, 
due to misalignments, thermal distortions of the beam and other error sources 
the computer generated commands will not be exactly correct. Hence, while 
in the autotrack mode the actual gimbal angles can be compared with the com- 
manded to produce error data. This data can then be fitted to seme selected 
model. See Figure 5.0-1. 



Figure $.0-1 Comparison of Actual and Commanded Gimbal 
Angles During Autotrack Mode 

The model which will be used in the simulation to represent the gimbal angle 
error is 


€ a = A a + B a sin % (t " t 0 ) + c a 008 % (t " t o ) 


= + B-, sin w (t-t ) + C fl cos a> (t-t ) 

p p p o o p o o 

where is the average orbit rate (-—) . 
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However, in order to illustrate the fitting procedure the following simple 
model will be used 


€ ■ + B., sin to (t-t ) 

a a a o o 


€ 0 " A 0 + B /3 8111 . W o (t-t o ) 


Defining x = 


A 


a 


B 


a 


(treating the axes independently) 



Vi>’ 


1 

i 

A 

Z "• 

w 

• 

• 

A 

s 

• 


€ (t ) 


• 

e 




_ a 
nr* 


1 sin to (t_-t ) 
oio 

1 sin to (t„-t ) 
o £ o 


1 sin to (t -t ) 
o n o 


x - H X 


The object is to find an £ ® % such that 

1 * a a A l A 

5<Z-HX, l-g*> » £<* , z> - J 


is a minimum. 


i.e. 


3j 


m A m 

* -a = (Z “ £ x) (- H) 

ax 


= cr 
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The above equation implies 

t <p A 

+f y s f HX 

Therefore if [H^ H] has rank two [H^ Hi" 1 exists and 

X - - Ch t H]' 1 H T jr 

It is suggested that during the actual calibration that the data vector £ be 
filled by taking points equally spaced around the orbit. 

The simulation will contain the capability of adding in the error model 
specified previously, but it will be up to the user to generate the coefficients 
for the best fit. The simulation can be used to generate the data from which 
the fit can be made once all the important error sources have been added to the 
program. Bias and thermal errors will be the error sources which will be mainly 
compensated. 

By doing this type of calibration it should be possible to reduce the overall 
pointing errors of the computer steering mode of operation. It will also 
reduce down times, i.e., if the monopulse momentarily looses track and the 
computer steering is immediately initiated, a reacquisition search mode require- 
ment will be less likely. 

Actual Dynamics, Kinematics, and Geometry Model 

This section describes the mathematical model used by the simulation program 
for the environment in which the antenna system is to operate. Portions have 
been extracted from the Task 2 Report — Mathematical Model for the Pointing 
Requirements of the Space Station High Gain Antenna - dated 28 Oct. 1972, 
and the Task 3 Report — Math Modeling of Influence Factors — dated 22 May 1972. 

The programming of this modeling is contained in subroutines KINMAT, ANTDYN, 
EPHEM1, EFHEM2, EIHEM3, SON, and M00N which are called by SETUP to determine 
the actual line of sight and sun vector directions. 
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3«U*2.0 Antenna Mounting Dynamics and Deflections (ANTDYN ) 

The antenna dynamics model has its origins in some work performed by 
G. Margulies . This model was modified for further work at LMSC. The model 
at that stage was present in Appendix D of the original proposal for the 
"Computer Simulation of Space Station Computer Steered High Gain Antenna" . 

The model contained within this section is a modification of those equations 
to meet the purpose and assumptions of this study. The resulting model is 
described in Section 3*1 of the Task 3 Report. 

In order to include the effects of the antenna boom flexibility (and space 
station flexibility in general), the dynamics equations could not be separated 
from the flexibility equations . Hence a great deal of effort was required to 
obtain the equations describing the overall rigid body and flexibility model. 
To accomplish this most efficiently, the aid of G. Margulies and his associate 
J. Ho (of the R and D division of LMSC) was enlisted. The work presented is 
a condensation of the work done by J. Ho in May of 1972. It is felt that this 
resulting model correctly models the configuration of the space station gim- 
balled antenna. However, there may be some difficulty in actually arriving 
at the required parameters associated with the flexibility model. 


G. Margulies, "Mathematical Model for the Attitude Dynamics of the Apollo 
Telescope Mount Satellite, LMSC /T -50-66-1, 12 December 1966. 
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3«U*2.1 Notations and Definitions 

Square brackets (C ]) are used here to indicate matrices with the symbol or 
* 

symbols inside labelling the whole matrix (one symbol) or the individual 
elements (appropriate number of symbols) or partitioned submatrices (appropri- 
ate number of symbols spearated by dashed lines). Vectors and diadics 
(operators) are indicated by symbols with cne or two (respectively) bars over 
them. When a vector or diadic symbol appears inside of a set of square 
brackets with a subscript P, or A this indicates the matrix of coordinates of 
the vector or diadic in the space station (P for platform) or antenna (a) 
reference frame. 


The symbols used here are defined in the following lists. 
Antenna Characteristics 



= inertia diadic of antenna 
= mass of antenna 

= displacement vector, antenna center of mass to 
hinge point 

^ unit vector in direction of antenna fixed gimbal axis 


Space Station Characteristics 



space station mass 

displacement vector, undeformed space station center of mass 
to hinge point 

unit vector in direction of space station fixed gimbal axis 

bending mode frequency for free vibration of non-spinning 
space station in mode described by p (below) 
bending deformation mode shape, a vector function of spatial 
coordinates which is analytic, has zero value and zero curl 
at the hinge point. 
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In the following definitions J ( ) dm indicates integration over the mass 

P 

distribution of the space station, and r^ is the displacement vector from the 
hinge point to the element dm . 

m,_ - /V dm 


v. = & 
h f* 


<x> » 


dm 


P P 


Y " ~ f (r x d 


m 


P P 

I - -J- 


ssjfr ■ A l 

rr - x 


~ j<p 2 E -<p <p | dm,E is unit diadic 


T A 
1 
C 


bending mode slope ^ y x ^ at control position 


Variables 


CO == angular velocity of space station at hinge point 

to = angular velocity of space station at control position 

c 

a. = gimbal angle about e a 

0 = gimbal angle about e« 

A 

q = bending coordinate in terms of which the translational and 
angular displacements of the space station due to bending 
are given respectively by q <p and q ^ V x tp 
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Forces 


T q = torque about gimbal axis e. 


= torque about gimbal axis e^ 


Derived Quantities 


R A 

«S 

L 


R 22 m B 


t Ce a 3 p ; Ceglp 3 

“t “ m A + m P 
®B = Ic f/ n 'T 

{ (r 2 . r 2 ) e - r 2 r 2 | 

*A + R 22 

_ A * - * - 

7? - a e a + P e p 

e = a P (e^x ep) + (a>x T 7 ) 

r[ = Rp - q<j> 

_ A * — _ -i 

b = 2 q <J> — CO x 

Q ^ — (co + 77 ) x L •■(<*>+ tj) — L . a - iigR^ x (Wx b) 

®12 = m B * R 2^ E ” **1 ** 2 } 

*=• A — 

N = L - 

$3 p = m B Ca 2 X <t> ] p 

„ A 2x2. 

M - m — m p <l> /m- 

A _ / ra t> \ 

C *V m B [Y + (m^ R 1 ~i^ R 2) X< * 3 P 
K * (m^-nip^ 2 ) w 2 


= — K q + 1 . W , (!P + q E ) .CO 


2 — — 

(nip/n^) <t> . (cox b) + m p <$ . j x a - Wx (wx 
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3.1*. 2. 2 Antenna Dynamics Equations 

For the special case in which the space station dynamics are not considered 
and the space station rate U ) , and acceleration u > are treated as known inputs 
the antenna dynamics equations may be expressed as 


r — 


— 

• • 
a 

[oO T [1] [0] i Cq] t ta 





Ji. 

Gfl T Co] ! m 


•• 

q 




[o] T Cq'] ♦ 


&]- 


Cg] t [S'] [«] 


$ - wl* ci] 


Eq. 3*0.1 


where the subscript P has been dropped from [G]p, [L]p, tQp, CQ 3 p > [N 3 p , 

U]p, and Cco]p both for clarity and to emphasize the fact that the equations 
are independent of the reference frame chosen for coordinatization. 

Equation 3*0.1 was obtained from Ref. 2 by straightforward translation to the 
notation of Ref. 1 followed by specialization to the antenna dynamics only 
case considered here. 


Since the space station is not rigid its angular velocity is not unique . The 

W which appears in Equation 3*0.1 is the space station angular velocity at the 

hinge point which is not in general the controlled part of the space station 

• • 

for which the angular velocity OJ and acceleration co are known. These 

c c 

quantities are related by characteristics of the bending mode shape: 


co 

c 


Z) + q T 


U> ** u> + q T 
c c 


Eq. 3-0.2 
Eq. 3*0.3 
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Since q is a state variable Eq. 3*0.2 may be used £ 0 compute a> from the input 
Zu . However q is not a state variable and is not available for the analogous 
computation of tf, hence Eq. 3*0.3 is combined with Eq. 3.0.1 to obtain 



fc] T Cl] to] ; to] T (cq - [S'] CT c 0) 


fee] 1 to3 ! ■ m - W T CT ] 

• C 


• * 

a 


• i 

0 



to] T K’ i * 


a 




- [g] t [n f 3 Ca> 


i- c*i T t«] 


Eq • 3-0. U 


These equations are used in the simulation program to compute the derivatives 

•••••• •• • * 

a , p, and q of the state variables a, P, and q respectively from the state 
variables and the driving quantities T a> T^, W , i,ud supplied by other 
elements in the program. 
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3«U»2.3 Bending Inputs for Sample Mode Shape 

As an example of a flexible S/S model consider the simple case of a rigid main 
body with a flexible antenna mounting boom depicted below. 




If the flexible boom is taken to be of negligible mass then the modal dis- 
placement <p is only important for the main body, where it may be expressed as 


(r) = 

C£L - 


A A A 

-aj-azj+ttik 


0 

-a -a z 
ai 


The bending inputs for this mode shape are evaluated as follows: 

\ - yv 2 dn = y'{ (a + az ) 2 + ( a *) 2 } 


dm = 


2 2 P 

s m p a + a - 
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a_ 


P a 

+ m “ 

F a 

4 






a 

o 

o 


for control system on main body. 


The mass distribution parameters (inertias) used in the above results are 
defined as: 




X dm, similarly for T & Z 


XI dm, similarly for XZ & YZ 


P A p p 


IJj - + I' , similarly for II & ZZ 


3*U«3 Actual Kinematics (K3NMAT) 


The mathematical model for obtaining the space station attitude from the space 
station angular rate and the orbital rate is identical to that used by the on 
board kinematics computations (Section 3>3«2). 
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3.1*.i* Vehicle Ephemerides 
Mathematical Discussion 

The equations used by subroutine EFHEM1 to compute the satellite ephemeris 
information are derived in the LMSC document The Motion of a Satellite About 
an Oblate Earth, H. W. Small, Nov. 10, 1961, LMSC/A086?5>6, and are more 
clearly presented in the LMSC document The Motion of a Satellite About an 
Oblate Earth — A Guide to the use of the Solution , P. Norov, Nov. 29, 1962, 
LMSC /A 31*021*6. The statement of these equations without the derivation is 
presented here. 

The position vector R and the inertial velocity vector V coordinatized in 
the earth reference frame are obtained from the instantaneous orbit elements 
P, A, B, i, ft, u by the relations 


£ § 

1 + A cos u + B sin u R 


V = + A cos u + ® si 11 uje^ + ^A sin u-B cos u^ e R J 

where jz = 1.1*0761*5 x 10^ ft^/sec^. 




cos u cos 0 — sin u cos i sin 0 

cos u sin 0 + sin u cos i cos O 

sin u sin i 

— sin u cos 0 — cos u cos i sin 0 

— sin u sin 0 + cos u cos i cos 0 

cos u sin i 


and. 
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The first five of these instantaneous orbit elements are obtained from the 
mean orbit elements Pg, eg, fcjg , ig, (which are constants) and u by the 
relations 

P - P B (1 * Z,^ sin 2 ijj) 

A = + B q G — J 2 C(1 - 6g) cos u (P 2 + sin 2 ig cos 2 u) + 

P« i (A + (1 + ♦ ) cos u] 

coo o 

B = B q — A q G — J 2 C(l - e^) sin u (P 2 — sin 2 ig sin 2 u) + 

P 9 ♦ (B + (1 + ♦ ) sin u] 

coo o 

1 " ig ‘ 5 J 2 slB ^ C 6s ij 


0 * (f J 2 cos ^ u + ^2 a 2 COS *3 


where 


j 2 4 j. 


4 i 


(^J = .0010827 (sssfmf 

2 2 3 3 

(cos u — sin u) + 2 (A cos J u — B siir u) 


A 3 . 

= 2 sin u 


3 3 

cos u + 3 B cos u — A sin u — B cos u 
o o o 


A^. = e B cos [ai B + \ J 0 (2 - 


.2 


B 2 2 


B 


eg sin 


in [w fi + | J 2 (2 - 


sin 


. 2 
sin 


V u;l 

ig) > 1 ] 


A -T 


G ■=’ Ojd — ^ sin^ ig) — ^ Jg (A sin u - B cos u) 


4 3 . 2 . .2 1 

= I sin ig sin u - 2 


I* = 1 + A coo u + B sin u 

0 0 o 
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while u is obtained by numerically integrating the equation 

u = ^1 + A cos u + B sin u) |l + 3 J 2 ("F") cos2 i 


sin u x 


c 


1 + A cos u + B sin 


in u^ 


The transformation T_ from the earth reference frame to the orbit reference 

o-I 

frame and the inertial angular velocity of the orbit reference frame cu 

are obtained from a 313 euler angle transformation with angles Q,, i, u and 
. . • & m 

rates fl, i, u where 

= n— <^,t, angular velocity of earth =.72921321 x 10 ^ sec 1 


u 


m 


i defined above 
u modulo 2 m 


* 

i 


= 3 '~2 GpO (l + A cos u + B sin uj c 

_ ,T> v2 P 1 


= u defined above 


1 + A cos u + B sin 


in u) 


:os i 


cos i 


. 2 
sin u 


cos u sin l 


followed by the interchange of axes given by the transformation 
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3. I*. 9 Sun Ephemeris 

Mathematical Discussion 


In the initialization section, Subroutine SUN first converts the input date 
(e.g., 062767 for June 27, 1967) to Modified Julian Date (MJD), which is 
defined as 

Modified Julian Date = Julian Day Number — 2,1*00,000.9, 


where the Julian Day Number is the number of days which have elapsed from 
Greenwich mean noon on January 1, 1*713 B.C. (Julian proleptic calendar), 
to Greenwich mean noon on the cale'ndar date in question. (Thus,, at noon 
June 27 , 1967, JD = 2,1*39,662, and MJD = 39*661. 9« ) The MJD is thus a mea- 
sure of the JD in which the "21*" is dropped and the day advances one unit at 
midnight rather than at noon. 

The mean orbit elements of the sun at midnight ( 0* 1 GMT) beginning the input 
date (or the following day if the seconds part of the base time is greater 
than 86,1*00) are computed from formulas referred to an epoch date of January 
1.0, 1900 (midnight beginning January 1, 1900,MJD 19020) rather than January 
0.9, 1900 (noon on December 31, 1899) as in Ref. 1, the Explanatory Supplement 
to the American Ephemeris and Nautical Almanac. The distinction between 
Greenwich Maan Time (or UT, Universal Time) and ephemeris time ET is ignored 
for these purposes. Let d be the number of ephemeris days which, have elapsed 
from the epoch date to the base midnight date: 

d = MJD - 19,020. (NDAY (3)) 


also, let T be this time period measured in Julian centuries of 36929 
ephemeris days: 

T = iSiZ (TCNTRY) 
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(TOWARD MEAN EQUINOX OF 
BASE MIDNIGHT DATE) 


Figure 6.1 The Inertial Coordinate System and Solar Angular Orbit Elements 
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Then the mean elements of the sun with respect to an inertial coordinate frame 
(Fig. l) referred to the mean equator and equinox of the base midnight date 
are given by the following expressions (cf. Ref. 1, p. 98): 

Mean anomaly Mg, measured from mean perigee of the sun's orbit, along the 
ecliptic to the mean sun: 

M 0 , B - ■ . 358 ? 96863013 + (1° - ? 011*399733) d 

- ? 0001$ T 2 - ? 000003 T 3 . (CM) 

Mean eccentricity e Q : 

♦ 

e„ _ = .0167510391* - .00001*180 T — .000000126 T 2 (ECC) 

Mean obliquity of the ecliptic €: 

= 23?1i$229382 - ? 013012$ T - ? 00000161* T 2 

+ ? 000000$03 T 3 . (0BL) 

Geometric mean longitude Lg, measured from the mean equinox of the base 
midnight date, along the ecliptic to the mean sun: 

L_ _ = 280? 189$0367 + (1° - ? 011*352661*6) d 

W yD 

+ ? 000303 T 2 , (CL) 

\ 

All of the angles are taken modulo 360° and converted to radians. 

In addition to the base solar elements, the Greenwich mean sidereal time 
(GMST), also referred to as the right ascension of Greenwich 0^ or 
(imprecisely) the Greenwich hour angle GHA is computed at the base midnight 
date: 

* 

u n _ = 100? 181181*91* + (1° - ? 011*352661*6) d. (CCMG) 

U ,D 

a;, is also taken modulo 36O 0 and converted to radians. 

U y b 
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On subsequent calls. Subroutine SUN computes the inertial and earth-referenced 
direction cosines of the true sun at any time T (in seconds) after the input 
base time of the run, or TIME = TIMERF +T seconds after the base midnight 
date: starting from the mean angular elements at this time, 

M_ = _ + 1.9909687222 x 10“ 7 

(!) (!) ,B 

L ft ■ L n „ + 1.9910638014 x 10" 7 

0 0,B 

the true longitude of the sun A^ is found by computing the equation of the 
center € from a series expansion in the mean anomaly (Ref. 2, p. 77) s 


TIME 

rad. 

(AM) 

TIME 

rad . , 

(AL) 


‘ l e + 4 e» 


rad. 


(ARGS) 


where 


= (2e — jj e^) sin M + (|j e 2 — e^) sin 2M 

,, , ir , , (DARG) 

+ si* 1 3M + — e 4 sin I 4 M 

r ,5 A 27 2 2L. 12£ 2 . 2 . 

= ^'2 + cos M ^ e son M) e cos M 


and e is written for e 


+ (2 + 3e 2 cos 2 M — ^ e 2 sin 2 M)] esiriM, 


0,B" 


h 5 

The error in truncating this series expansion at the e term is of order e , 
-9 

or 1 x 10 radians. The inertial direction cosines of the true sun are 
then 


(SI) 


I 

M 

© 

, ^ 


cos A 0 

m_ _ 


sin A cos €L 

0,1 


0 B 

1 

M 

*\ 

© 

C 1 


cos Ag sin 
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to-inertial transformation matrix is used: At TIME seconds after the 


To find the direction cosines in the earth-fixed coordinate system, the earth 
to-inertial transformation matrj 
base midnight date, the GMST is 

°B = °b,B t ‘k’ TIME 


rad . , 


(OMj) 


where is the inertial rotation rate of the earth, 

= 7 . 292 ll^Iii 7 x 10~* rad/sec. 

Neglecting precession and nutation effects, Tgj may be written 

cos -sin 0 

sin . cos 


l EI 


0 

1 


(OMjE) 


(ETI) 


so that the earth-referenced direction cosines of the true sun are 


O 

* 

1 


~^0,I ~ 

T 


m 0,E 

. = (t ei } 

*"0,1 

n 0,E 




(SE) 


The magnitude of the radius vector from the earth to the sun is then cal- 
culated. In astronomical units, it is given by 


, = (1-e*) 

l 0 1 + e cos (M 0 + . 


(PAU) 


Subroutine SUN also determines whether the satellite is in sunlight or not. 
The conditions which must be satisfied for the satellite to lie in the 
(cylindrical) shadow of the earth are (Fig. 2) 
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Figure 6.2 Geometry for Satellite in Earth's Shadow 



Figure 6.3 Orbit Beta Angle 
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cos < 0 


R sin t < 


^EQ 


where RgQ is the earth’s equatorial radius, R is the satellite radial distance, 
and t is the geocentric angle between the sun and the satellite. If u Q is the 
earth-frame unit vector in the solar direction, then the above conditions can 
be written in the vector form 


R . u. < 0 
0 


Rxu 0l < 


The integer NIGHT is set to 1 if the satellite is in shadow and 0 if it is 
in sunlight. 


Finally, SUN calculates the "orbit beta angle" fi (Fig. 3), which is the 


angle between the solar direction unit vector u @ and its projection on the 
satellite orbital plane. If C is the angle between u 0 and a vector normal 
to the orbit plane, 

sin £ = |£ 0 x v^j , (SBETA) 


then 


fi = 90° - C • 


(BETA) 
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' 3.U*6 Moon Ephemeris 
Mathematical Discussion 

In the Initialization Section, Subroutine MOON computes the mean orbit 

elements of the moon at the base midnight date (see SUN Mathematics Section). 

If T is the time elapsed from an epoch date of January 0.5, 1900 (noon on 

December 31, 1899), 

T - , (TCNTRY) 

then the mean elements with respect to the inertial coordinate system 
described earlier (Fig. 1 of Subroutine SUN section) are given by the follow- 
ing expressions (Ref.: Explanatory Supplement to the American Ephemeris and 

Nautical Almanac): 

Mean anomaly M^ , measured from the mean lunar perigee along the orbit to 
the mean moon: M^ =<f— T , 

M tf B = 296!l0U608 + 1*77,198°8U9108 T 

+ !o09190 T 2 + ! 0000032 T 3 . (AMI) 

Longitude of the mean ascending node , measured from the mean equinox 
of date along the ecliptic : 

fi <r,B = 2 ^* l8 3 2 75 - 193l|!lli201 T 

+ ! 002077 T 2 + .000002 T 3 . (iZMl) 

Mean argument of perigee measured from the mean ascending node along the 

t 

lunar orbit: a) 4 = T — , 

B = 75°lli628l + 6003°1760Ul T 

- !oi2li02 T 2 - !0000l5 T 3 . (WMl) 
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These angles are all taken modulo 360° and converted to radians. In addition, 
the following mean elements are considered as constants: 

Mean eccentricity e^ 


e <r,B 

. 051:9001*89 

(ECCM) 

Mean inclination i<j : 



sin B 

= .089683UUS 

(S3NI) 

COS 

= .99^97032 

(COSI) 

2 

Mean semi-latus rectum, = aj (1 - e^ ): 


p a,B 

1.257355 x 10 9 ft. 

(M) 


On subsequent calls. Subroutine MOON computes the distance and inertial 
coordinate frame conponents of the true moon at any time T (in seconds) 
after the input base time of the run, or TIME = TIMERF + T seconds after 
the base midnight date; starting from the mean angular elements at this time, 


M < = M «B + * 228o2 7135 • T rad. 

(AM) 

B -.000921:22029. t rad. 

(0M) 

03* = Wj B +.0028685883 • Trad. , 


t = TIME/86, 1:00, 

(TIME) 


The true anomaly of the moon 6^ is found by computing the equation of the 
center from a series expansion in the mean anomaly (see SUN Mathematics 
Section): 

* 

e <L = M <j + ^ , rad (THM) 
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where 


’ 2e «, B Sln M « + t e «,B si " 2 «« • 

2 

The error in truncating this series expansion at the e term is of order 
3 

e , or .0002 radians. 


The magnitude of the geocentric radius vector is then 


R a 1 + e 


SiL 


«,b cos e c ' 


(RMS) 


and its components in the inertial coordinate system are 


L «,I 






where 


cos cos u ^ — cos i^ g sin 0^ sin u^ 

(sin ft* cos u * + cos i* _ cos ft* sin u*s ' c 

a 4 C,B ( <t) cos fc 


B 


— sin i - _ sin u* . c 
C,B <1 sin t 


B 


(sin ft^cos u^j + cos i^ B cos ft^ sin u^) sin 


+ sin i^ B sin u^ cos 


u « ’ e t* u s 


(UM) 


and is the mean obliquity of the ecliptic., calculated in Subroutine SUN. 
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APESNDIX A 

DECK SETUP AND INPUT LISTS 



This appendix describes the format and ordering of the data cards which are 
read by the program to enter the input data. 

Data cards are placed in the deck immediately after the program execution 
card (VXQT L00P for the complete simulation program or V XQT F0INT for 
the pointing requirements program) and are arranged as follows: 



VXQT P0INT 
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If more than one case is to be run, one or more blocks of the following might 
be inserted between the blank card and the "J" card (at *). 

/ ( blaW WBT- ~~~ 

/ PLOY cWangK <! CANNOT DCUIt Plots, ONLY MODIFY > 

/ L A i P { l: ChAt‘-.(*E. h kill CASE! 1 

^TITIX "CAST ” ■ 


r. 


000 0 00000 00000000000000000000000 OOOQOOOOOOOOOOOOOOOOOOOOOO 000 

miii mnniiniingii util ;o ?i n n h » » v a n » ji u « » n * u ■ « « »i <! «j « « « « a « vmi v u » « v m « «> u « 

1 1 1 n m n> 1 1 1 m n 1 1 1 1 1 1 1 1 1 1 1 m 1 1 1 1 1 m 1 1 1 1 1 1 m 1 1 1 1 1 1 1 m i m 

22222222222. 2222222222222222222222222222222222222222222222222222 


3 3 3t • 3l 3 3| 3 3 3 3 3 3 3 3 3 3 3 3 3 3 i 3 3 3 3 3 J 3 3 3 3 3 3 3 3 3 3 3 3 J 3 3 J 3 J 3 J 3 3 3 3 3 3 3 3 3 J 3 3 3 

444444444444444444444444444444444444(4444444(4444444444444444444 


To change the "nominal” data, i.e. clear V -ARRAY and FHEAD information before 
in put in g the next set, place a ”2 " card, (a 2 in column one) before the desired 
title card. 


CARD FORMATS 
Title 

Column 1 must be blank (see General Comments at end for extended feature) 
Data for V-array 

The data for the V-array is input in free field format: 



where 1 i is a V-array position index and is the data value associated 
with V(l^) . Blank(s) (at least one) must separate the 1^ and v^, and 1^ 
need not start in column one. 

The 1 are in integer format with no included blank(s) allowed since blank(s) 
are separators. 


A-2 


The v i may use any of the following formats: 

Integer - no decimal point in number 
Fixed point - value has a decimal point 

Exponential - number followed by the exponent as in floating point number. 
All the following are valid and equivalent. 

U2.5, U-2$ +1, U.25E1, U25E-1 

Example: 

! ,’!•-> J li 4 Pi 1 «r. . 1 


00000000000 8000000000 000 00000000000000000000000000000000000000000000000000000 

i i i 4 i ( i a * it m » u ii u it u ti n it n » » n » x n it n a ii v u m k » u a » n> n « o m « it n a it m it ;i m m » u n u » u ti ti u ii ci u ca n ii n u ii n it ii n it a 

ii iiiiiiiiinii 1 1 1 1 1 1 1 1 inn ini it 1 1 1 1 1 1 1 1 1 1 n 1 1 1 i i ii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

22222222222 222222222222222222 2222222222222222222222222222222222222222222222222 

33333 33333333 3333 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

44444444444444444444444444444 44444444444444444444444444444444444444444444444444 
5 5 5 S 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 

660666 66C6666G 666666 6666666 666666666 6666666666666 666 6S666666666666666666666666 

I 7 7 7 7 7 7 7 7 7 7 7 7 7 I 7 7 / 7 7 7 7 7 7 7 7 7 ? 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 
| 888888 88888888 8888888888688888888888 88888888888888888888888888888888888888888 
99939099999999999993909939999999999999999999999999999999999999999999999999999999 

I 7 3 4 . 4. I • * to III/ 13 ‘4 lb It 1/ ll 19 TO 21 Z7 74 ft ft 21 a ft a II 17 ’S » * / V a B « 41 47 «3 41 4S 4| 4} 43 X SI U X Ft ./ V X It t? tl fi# U ft 1/ X M ft J1 77 U 74 ft ft 77 ft * B 
flH VWl 

**** DATA MUST BE PAIRED, i.e., EVER! 1 ± MUST HAVE A v ± **** 

PLOT DATA 

Plots Iran the DEADAM package are available in four forms (see Appendix A). 

The first three are closely related since their formats are derived from the 
same array (FREAD). For discussion, they will be referred to as Full frame, 

1/3 frame , both produced on the FR-80, and printer plots plotted on the high 
speed printer . 

i 


At3 



These three plot forms are specified in the following manner. 
Full and 1/3 frame - 



P 1 d l p 2 d 2 * * * P i d i 


; i 


1 , 10 



is the plot number 

is the channel information for plotting on plot p^ 
cL may be any one of the following 

OOY or Y will plot output channel Y vs time 
X00 will plot time vs output channel X 

XY will plot X vs Y 

X and Y must be 2 digits (i.e., 21, 02, etc.) 


To designate that the first group of plots requested are to be 1/3 frame 
(three plots per frame) use 


f P 1000 n 


where n is the number of plots to be plotted three to a frame. All plots 
numbered greater than n will be full frame. If n is not an integral multiple 
of 3, there will be blank 1/3 frame plots an the ([j] + l)^ 1 frame. 


Printer plots are rough plots received as part of DEADAM printed output. To 
get printer plots in addition to FR80 plots, use the following data card: 



where m is a plot number specified on a "P" card that is to be print-plotted 
the n^ need not be sequential. 


k-k 




Plot Scale Card 


Those cards are again associated with the {{Lots specified on a "F 1 card and 
will therefore affect both FR80 and printer plots. They allow scaling to 
specified limits for these plots. 

Format 


S n 0 1^ . . . 1 A 

Where S specifies Plot scale card 

n number of plot to be scaled (corresponds to a "P" card 

specification) 

0 dummy (not used) 




1 . 

x 


limit being set (must be grouped: if X min set, X max 

must be set) according to following table. 



3 

h 

5 

6 


limit 

X min 
X max 

Y min 

Y max 


Data generated outside of specified limits will not 
be plotted . 


v. 

x 


limit value. Fairs (l.„ v ) have the same format as 
(V-array position, data value) pairs. 


Multiplot card 

Multiplot allows the plotting of up to 10 channels vs the sane abscissa; the 
plot is full frame. The multiplot card has the form 




n k 1 



\ + 1 


A -5 



specifies multiplot card 
multiplot # 

number of ordinate channels 

channels to be plotted on ordinate, 1 = 1, . . ., k 
channel to be used as abscissa 

(If 1^ + ^ is 0 or blank, l^ +1 is assumed to be time.) 

Scaling is automatically set. 

General Comments : 

A (@) punch on a Data or Plot cards will cause the program to ignore those 
columns following the (@) punch cn that card. At least one set of data informa- 
tion must have been set before the (@) punch. 

Title card extension ; If there is an in column 80, that card will be listed 
immediately; and the program will read the next card to get the heading informa- 
tion. This will continue until the last Title card read has no in column 80. 
This allows the user to append seme comments before his V -ARRAY list out. 

Entries in the V-ARRAY (entered via the above described data cards) are interpreted 
by the individual subroutines of the program according to the following input lists. 
The absolute location of an input variable is its V-ARRAY index. Relative loca- 
tions are also given since the starting index far the input data of each routine 
(except ADAMS and ANTIN) may be changed (see input list far subroutine SETUP). 


where M 
n 
k 


1 . 

x 


k + 1 
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Input List, For Subroutine SETUP/L00P and SETUP/P0INT 

Variable Relative Absolute Nominal 


Name 


Location Location 

Description 

Units 

Value 

NSS0RB 


1 

input index for EPHEM3 

— 

50 

NR0RB1 


2 

input index for EHEM1 

— ■ 

60 

NR0RB2 


- 3 

input index for EHiEM2 

— 

70 

NVSUN 


h 

input index for SUN 

— 

80 

NVSS1 


5 

input index for 3NTSS1 

- 

90 

NADYN 

* 

6 

input index for ANTBYN 

— 

600 

NCNTL7 

* 

- 7 

input index for CNTRL7 

— 

i 1*1 

NSERV0 

* 

8 

input index for SERV0 

— 

181 

NMNPLS 

* 

9 

input index for M0NFLS 

— 

201 

N0BR10 

* 

• - 10 

input index for 0BEPH1 

— 

560 

N0BR20 

* 

11 

input index for 0BEPH2 

' — 

580 

N0BSS0 


12 

input index for 0BEPH3 

- 

Soo 

NTDEF 

* 

13 

input index for THERMD 

- 

250 

NKNMT 

•5S- 

iU 

input index for KINMAT 

- 

300 

N0BKMT 

* 

15 

input index for 0BKNMT 

— 

325 

NSSRTE 

* 

16 

input index for SSRATE 

- 

330 

NCLBRT 

* 

17 

input index for CLBRTN 

— 

31*0 

GSLAT 

*«• 

21 

ground station latitude 

deg 


GSL0N 


22 

ground station longitude 

deg 


GSALT 

** 

23 

ground station altitude 

N.M. 


KHIINT 


Ul 

output option index (not used) 

— 

1 

KDEBUG 


1*2 

/ o -♦debug outputs each page 

— 

0 

KRELAY 


U3 

relay to be tracked 


1 

K0PT2 


1*1* 

option index (not used) 

— 

— . 

K0PT3 


1*5 

option index (not used) 

— 

— 

* — ► SETUP/L00P only 




-** — - SETTJP/P0INT only 
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Input List For Subroutine ADAMS (internal integration routine) 


Variable 

Relative 


Absolute 




Nominal 

Name 

Location 


Location 

Description 

Units 

Value 


INPADM 



1*00 

not used 


— 

- 

I0IN 

IN PALM 

+ 

1 

1*01 

[discontinuities at 

— 

0 






TDISC(l) 




I0MAX 

INPADM 

+ 

2 

1*02 1 

|l = IJZJIN, ... I0MAX 

— 

0 

TIGN 

INPADM 

+ 

3 

1*03 

initial on time 

periodic 

' — ' 

0 







time 



TINC 

INPADM 

+ 

h 

1*01* 

period 

1 discon- 

. — 

0 

TC0MAN 

INPADM 

+ 

5 

U05 

on duration 

tinuities 



IGUESS 

INPADM 

+ 

6 

1*06 

- 1 to suppress output at 

- 

0 


1 




disc . 





INPADM 

+ 

7 

1+07 

> 0 to suppress output at 

sec . • 

0 






disc, except for V(l*06) 








<T< V(l*07) 




TL 

INPADM 

+ 

8 

1*08 

initial value of 

time 

sec . 

0 

TINT 

INPADM 

+ 

9 

1*09 

output interval 


sec . 

1 

TMAX 

INPADM 

+ 

10 

1*10 

maximum value of 

time 

sec. 

i£o 

NIMAX 

INPADM 

+ 

11 

i*n 

maximum # integration 

— 

131071 * 






steps 




I0UTM 

INPADM 

+ 

12 

1*12 

maximum # output lines 

- 

2-9* 

DTIN 

INPADM 

+ 

13 

U13 

initial integration step 

sec. 






size 




DTMIN 

INPADM 

+ 

Ik 

UHi 

minimum integration step 

sec . 

2 *9 * 






size 




DTMAX 

INPADM 

+ 

15 

Ul5 

maximum integration step 

sec . 

2+5 * 






size 



. 23 V 6 * 
■1,6V 5 * 

ADLL 

INPADM 

+ 

16 

1*16 

terror test limits 

- 

ADUL 

INPADM 

+ 

17 

U17 

1 


— 

EFC0N 

INPADM 

+ 

18 

1*18 

predictor error limit 

— 

1 CT S » 

IDEBIN 

INPADM 

+ 

19 

U19 j 

1 print debug outputs for 

- 

0 

IDEBF 

INPADM 

+ 

20 

1*20 ! 

I IDEBIN < NI < IDEBF 

- 

0 


* values set internally if input is 0 
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Input List For Subroutine ADAMS (internal integration routine) Continued 


Variable 

Relative 

Absolute 



Nominal 

Name 

Location 

Location 

Description Units 

Value 

3ERR 

IN PALM + 21 

1*21 

1 to make predictor error 
test 

— 

0 

L0UTA 

INPALM + 22 

1*22 

1 to suppress plot output 
at disc. 

— 

0 

IITMAX 

INPADM + 23 

1*23 

# variable output incre- 
ments 

— 

0 

NXL 

INPALM + 2 h 

1*21* 

initial output channel 

— 

1 

NXH 

INPALM + 2 5 

1*25 

final output channel 

— 

21* 

TI0(l) 

INPADM + 26 

U26 

1st interval output spac- 
ing 

2nd interval output spac- 

sec. 

— 

TI0(2) 
TI0( 3) 

INPADM + 27 

1*27 

sec. 

— 

ing 

INPADM + 28 

1*28 

3rd interval output spac- 
ing 

sec. 


TI0(U) 

INPADM + 29 

U29 

l*th interval output spac- 
ing 

sec. 


TI0(5) 

INPADM + 30 

1*30 

5th interval output spac- 
ing 

sec . 



T0U(l) 

INPADM + 31 

1*31 

max. time for 1st interval 

sec. 

— 

T0U(2) 

INPADM + 32 

1*32 

max. time for 2nd interval 

sec. 

- 

T0U(3) 

INPADM + 33 

1*33 

max. time for 3rd interval sec. 

— 

T0U(U) 

INPADM + 3k 

1*31* 

max. time for l*th interval 

sec. 

- 

T0U(5) 

INPADM + 3$ 

1*35 

max. time for 5th interval 

sec . 

— 

TDISC(l) 

• 

INPADM + 36 
• 

1*36 

• 

discontinuity time 

• 

sec. 

• 


• 

• 

TDISC(29) 

• 

• 

INPADM + 6U 

• 

1*61* 

• 

• 

discontinuity time 

• 

• 

sec. 

— 

ADDCm(l) 

INPADM + 6 5 

1*65 

additive constant 1 

— 

10 

ADDC0N(2) 

INPALM + 66 

1*66 

additive constant 2 

— 

10 

ADDC0N(3) 

INPADM + 67 

1*67 

additive constant 3 

- 

10 

ADDC0N(U) 

INPADM + 68 

1*68 

additive constant 1* 

— 

1* 

ADDC|ZN(5) 

INPADM + 69 

1*69 

additive constant 5 

- 

1* 
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Input List For Subroutine ADAMS (Internal integration routine) Continued 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 


Units 

Nomina 

Value 

addc#n(6) 

IN PALM + 70 

1*70 

additive 

constant 

6 

— 

I* 

ADDCtfN(7) 

INPADM + 71 

1*71 

additive 

constant 

7 

— 

l* 

ADDCI2(N(8) 

3NPADM + 72 

1*72 

additive 

constant 

8 

— 

1* 

ADDC0N(9) 

INPADM + 73 

1*73 

additive 

constant 

9 

— 

1* 

ADDCdN(lO) 

INPADM + 71* 

1*71* 

additive 

constant 

10 

— 

1* 

ADDC#N(ll) 

INPADM + 75 

1*75 

additive 

constant 

11 

— 

1* 

ADDCdN(l2) 

INPAEM + 76 

1*76 

additive 

constant 

12 

— 

.5 

ADDC2lN(l3) 

INPADM + 77 

1*77 

additive 

constant 

13 

— 

.5 

ADDC0N(ll*) 

INPADM + 78 

1*78 

additive 

constant 

H* 

— 

10 

ADDCgjN(l^) 

INPADM +79 

1*79 

additive 

constant 

15 

— 

l 

ADDClZN(l6) 

INPADM + 80 

1*80 

additive 

constant 

16 

- 

10 

ADDCJZ!N(17) 

INPADM + 8l 

1*81 

additive 

constant 

17 

- 

1 

ADDC2lN(l8) 

INPADM + 82 

1*82 

additive 

constant 

18 

. — 

.005 

ADDCgftl(l9) 

INPADM + 83 

1*83 

additive 

constant 

19 

- 

.005 
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Input List For Subroutine ANTIN 




Variable Relative Absolute 

Name Location Location 

Description 

Units 

Nominal 

Value 


FLAG 

mxv 


no 

flag: if 0 ANTIN not 
caned 

mm 

0 

BMW 

EDXV 

+ l 

111 

beam width factor 

deg 

70 

CSN 

hkv 

+ 2 

112 

feed option parameter 

— 

2.5 

DBK 

IDXV 

+ 3 

113 

dia. of aperture blockage 

in 

0 

DIA 

IDXV 

+ k 

11U 

reflector diameter 

in 

60 

EFF 

IDXV 

+ 5 

n5 

aperture efficiency 

% 

55 

FTV 

IDXV 

+ 6 

116 

distance feed to vertex 

in 

0 

FVD 

I3?XV 

+ 7 

117 



•U 

GFE 

IDXV 

+ 8 

118 

peak gain of feed 

— 

0 

GHZ 

IDXV 

+ 9 

119 

operating frequency 

G hz 


HUB 

IDXV 

+ 10 

120 

radius of reflector "hub" 

in 

0 

M0N 

IDXV 

+ 11 

121 

not used 

- 

— 

PAN 

IDXV 

+ 12 

122 

# of reflector panels 

- 

0 

PFW 

IDXV 

+ 13 

123 

integration step(pts/ 
wavelength) 

— 

2 

SIM 

IDXV 

+ 1U 

12U 

flag: 0 for antenna study 
only 

— 

1 

SMP 

IDXV 

+ 15 

125 

flag: 1 for simple model 

— 

0 

XAP 

IDXV 

+ 16 

126 

\ 


— 

XAT 

mxv 

+ 17 

127 

/ 


_ 

XVP 

mxv 

+ 18 

128 

)not used 


; — 

XVT 

mxv 

+ 19 

129 

) 


— 

P0L 

mxv 

+ 20 

130 

polarization flag: 1 for 


1 


linear dipole, 2 for 
linear slot mode, 3 for 
right hand circular, and 
U for left hand circular 
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Input List For Subroutine ANT3N Continued 


Variable Relative Absolute 

Name Location Location Description Units 


me 

IDXV 

+ 21 

131 

For SIM = 0 
constant (p \ 


deg 

THS 

3DXV 

C\J 

CO 

+ 

132 

initial 8 / 

! 

•for 8 scan 

deg 

THD 

IDXV 

+ 23 

133 

8 increment ( 

deg 

THF 

IDXV 

+ 21* 

13U 

final 6 : 

1 

deg 

THC 

IDXV 

+ 25 

135 

constant 8 


deg 

ms 

IDXV 

+ 26 

136 

initial <p J 

I 

►for <p scan 

i 

deg 

PHD 

IDXV 

+ 27 

137 

<p increment | 

deg 

PHF 

IDXV 

+ 28 

138 

final (p j 

1 

deg 


Nominal 

Value 
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Input List For Subroutine M0NPLS 


Variable 

Relative 


Absolute 



Nominal 

Name 

Location 


Location 

Description 

Units 

Value 

M0N 

INDEXV 


201 

option index: 0-*- no 
monopulse error signals, 
it — ► four horn system, 

$ five horn system 


It 

GDB 

INDEXV + 

1 

202 

peak antenna gain 

db 

0* 

P 

INDEXV + 

2 

203 

[beam shape parameters 

db/deg^ 

3-0* 

T 

INDEXV + 

3 

2 Oil 

\ 

db/deg 

3.0* 

US 

INDEXV + 

h 

20$ 

horn spacing 

deg 

.5 

VS 

INDEXV + 

$ 

206 

horn spacing 

deg 

.5 

AS 

INDEXV + 

6 

207 

monopulse gain factor 

— 

.16 

THRESH 

INDEXV + 

7 

208 

monopulse threshold gain 


.5 

EL 

INDEXV + 

8 

209 

[error signal gains 

— 

.16 

K2 

INDEXV + 

9 

210 

) 

— 

.16 

THETAS 

INDEXV + 10 

211 

monopulse activation cone 

deg 

1.0 


angle 


indicates values that may be computed by ANTIN 
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Input List For Subroutine THERMO 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 

K0FT 

IDXV 


2$0 

option index (not used) 

- 

— 

DIAM 

mxv + 

1 

251 

antenna diameter 

ft 

30. 

FREQ 

hjxv + 

2 

252 

antenna operating fre- 
quency 

G Hz 

8.25 

AO 

IDXV + 

3 

253 


— 

1 

A1 

IDXV + 

a 

25U 


— 

0 

A2 

IDXV + 

5 

255 


' — 

0 

BO 

IDXV + 

6 

256 


. — 

1 

B1 

IDXV + 

7 

257 


■ — 

0 

B2 

IDXV + 

8 

258 


— 

0 

CO 

IDXV + 

9 

259 


— 

1 

Cl 

IDXV + 

10 

260 


— 

0 

C2 

IDXV + 

11 

261 


— 

0 

DO 

IDXV + 

12 

262 


■ - 

1 

D1 

IDXV + 

13 

263 


— 

0 

D2 

IDXV + 

1U 

26^ 


— 

0 

EO 

IDXV + 

15 

265 


— . 

1 

El 

IDXV + 

16 

266 


- 

0 

E2 

IDXV + 

17 

267 


— 

0 

GO 

IDXV + 

18 

268 


— 

1 

G1 

IDXV + 

19 

269 


— 

0 

G2 

IDXV + 

20 

270 


— 

0 

CDHII 

IDXV + 

21 

271 

) constant (indep. of sun 

deg 

0* 

CDTHA 

* 

IDXV + 

22 

272 

(vector) beam shift angles 

deg 

0* 


indicates values that may be computed by ANTIN 




Input List For Subroutine SSRATE 

Variable Relative Absolute Nominal 


Name 

Location 

Location 

Description 

Units 

Value 

K0PT 

IDXV 

330 

option index 

- 

0 




For K0FI = 0 S/S rate is 
set equal to orbit rate, 
no inputs 






For K0PT - 1 



gtopip(i) 

IDXV + 1 

331 

S/S X-axis rate (constant) deg/sec 

— 

0MPIP( 2 ) 

IDXV + 2 

332 

S/S Y-axis rate (constant) deg/sec 

— 

0MPIP(3) 

IDXV + 3 

333 

S/S Z-axis rate (constant) 

deg/sec 

— 




For K0PT = 2 



pg&ipip(i) mxv + i 

331 

S/S X-axis accel. 
(constant) 

deg/sec 2 


P0MPIP(2) IDXV + 2 

332 

S/S Y-axis accel. 
(constant) 

deg/sec 2 

— 

P0MPIP(3) IDXV + 3 

333 

S/S Z-axis accel. 
(constant) 

deg/sec 2 


R1 * 

IDXV + U 

33U 

S/S X-axis rate 
(initial) 

deg/sec 2 


R2 

IDXV + $ 

335 

S/S Y-axis rate 
(initial) 

2 

deg/sec 

— 

R3 

IDXV + 6 

336 

S/S Z-axis rate 
( initial ) 

deg/sec 2 

— 
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Input List For Subroutine KINMAT 


Variable 

Name 

Relative 

Location 


Absolute 

Location 

Description 

Units 

Nominal 

Value 

IR0T 

INDEXV 



300 

input rotation sequence 

■ — 

321 

FHI(l) 

INDEXV 

+ 

1 

301 

S/S X-axis rotation angle 

deg 

- 10 

HII(2) 

INDEXV 

+ 

2 

302 

S/S Y-axis rotation angle 

deg 

0 

HU(3) 

INDEXV 

+ 

3 

30 3 

S/S Z-axis rotation angle 

deg 

0 

N0ER(2tr 

INDEXV 

+ 

19 

319 

# of optional output 
sequences 


0 

RCHNG 

INDEXV 

« 

+ 

20 

320 

optional output rotation 
sequence 



RCHNG 

INDEXV 

+ 

21 

321 

optional output rotation 
sequence 

■“ * 


RCHNG 

INDEXV 

+ 

22 

322 

optional output rotation 

- 

— 


sequence 
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Input List For Subroutine EHIEM 3 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

De scripts' on 

Units 

Nominal 

Value 

K0PT 

INDEXV 

50 

input option index 

— 

3 




For K0PT = 1 



GMEGA0 

INDEXV + 1 

51 

average orbit rate 

rad/sec 

— 

EB 

INDEXV + 2 

52 

orbit eccentricity 


— 

CLINB 

INDEXV + 3 

53 

orbit inclination 

deg 

— 

WB 

INDEXV + k 

5U 

perigee argument 

deg 

— ■ 

AMU 

INDEXV + 5 
1 

55 

initial true anomaly 

deg 

— 

CAPWB 

INDEXV +6 

56 

ascending node longitude 

deg 

— 




For K0PT =2 



IERIOD 

INDEXV + 1 

5i 

nodal period 

sec 

— 

IRAD 

INDEXV + 2 

52 

perigee radius 

nm 

— 

CLINB 

INDEXV + 3 

53 

orbit inclination 

deg 

— 

WB 

INDEXV + U 

5U 

perigee argument 

deg 

— 

U 

INDEXV + 5 

55 

initial satellite argument 

deg 

■ - 

CABJB 

INDEXV + 6 

56 

ascending node longitude 

deg 

— 




For K0FE =3 



R0MNAK 

INDEXV + 1 

51 

orbit semi -major axis 

nm 

381*3.1 

EB 

INDEXV + 2 

52 

orbit eccentricity 

— 

0. 

CLINB 

INDEXV + 3 

53 

orbit inclination 

deg 

90. 

WB 

INDEXV + U 

$h 

perigee argument 

deg 

90. 

AMU 

INDEXV + 5 

55 

initial true anomaly 

deg 

- 90. 

CAPWB 

INDEXV + 6 

56 

ascending node longitude 

deg 
• - 

90. 
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Input List For Subroutine ETHEM 1 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

V 

Description 

Units 

Nominal 

Value 

K0PT 

INDEXV 

60 

input option index 

- 

3 




For K0FT = 1 



OMEGA 0 

INDEXV + 1 

6l 

average orbit rate 

rad/sec 

— 

EB 

INDEXV + 2 

62 

orbit eccentricity- 

- 

— ' 

CLINB 

INDEXV + 3 

63 

orbit inclination 

deg 

— 

m 

INDEXV + h 

6U 

perigee argument 

deg 

— 

AMU 

INDEXV + 5 

65 

initial true anomaly- 

deg 

— 

CAFWB 

INDEXV + 6 

66 

ascending node longitude 

deg 

— 




For K0FT = 2 



PERIOD 

INDEXV + 1 

61 

nodal period 

sec 

— 

RIAD 

INDEXV + 2 

62 

perigee radius 

nm 

— 

CLINB 

INDEXV + 3 

63 

orbit inclination 

deg 

■ 

V© 

INDEXV + U 

6U 

perigee argument 

deg 

— 

U 

INDEXV + 5 

65 

initial satellite argument 

deg 

— 

CAFWB 

INDEXV + 6 

66 

ascending node longitude 

deg 

- 




For K0FT =3 



R0MNAK 

INDEXV + 1 

61 

orbit semi -major axis 

nm 

22??0.8 

EB 

INDEXV + 2 

62 

orbit eccentricity 

- 

.00011 

CLINB 

INDEXV + 3 

63 

orbit inclination 

deg 

.01 

WB 

INDEXV + h 

Sh 

perigee argument 

deg 

0 

AMU 

INDEXV + 5 

65 

initial true anomaly 

deg 

0 

CAR© 

INDEXV + 6 

66 

ascending node longitude 

deg 

15. 
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Input List For Subroutine EFHEM 2 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 

K0PT 

INDEXV 

70 

input option index 

— 

3; 




For K0PT - 1 



OMEGA0 

INDEXV + 1 

71 

average orbit rate 

rad/sec 

— 

EB 

INDEXV + 2 

72 

orbit eccentricity 

— ■ 

— ' 

CLINB 

INDEXV + 3 

73 

orbit inclination 

deg 

— 

V® 

INDEXV + U 

74 

perigee argument 

deg 

— • 

AMU 

INDEXV + 5 

75 

initial true anomaly 

deg 

— 

CAR© 

INDEXV + 6 

76 

ascending node longitude 

deg • 

— 




For K0FT - 2 



PERIOD 

INDEXV + 1 

71 

nodal period 

sec 

— 

FRAD 

INDEX V + 2 

72 

perigee radius 

nm 

— 

CLINB 

INDEXV +3 

73 

orbit inclination 

deg 

— 

WB 

INDEXV + 4 

74 

perigee argument 

deg 

— 

U 

INDEXV + 5 

75 

initial satellite argument 

deg 

— 

CAR© 

INDEXV + 6 

76 

ascending node longitude 

deg 

— 




For K0FT = 3 



P.0MNAK 

INDEXV + 1 

71 

orbit semi -major axis 

nm 

22770.8 

EB 

INDEXV + 2 

72 

orbit eccentricity 

— 

.00011 

CLINB 

INDEXV + 3 

73 

orbit inclination 

deg 

.01 

V/B 

INDEXV + 1* 

74 

perigee argument 

deg 

0 

AMU 

INDEXV + 5 

75 

initial true anomaly 

deg 

0 

CAR© 

INDEXV + 6 

76 

ascending node longitude 

deg 

145. 
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Input List For Subroutine SUN 



Variable 

Name 

Relative 

Location 

Absolute 

Location Description 

Units 

Nominal 

Value 

DATE 

TM0D 

INDEXV 
INDEXV + 1 

80 date (month day year) 

81 initial time of day 
(Greenwich) 

sec 

032172 

21600 


Input List For Subroutine INTSSI 


Variable Relative Absolute Nominal 

Name Location Location Units Value 


A 

3DXV 



90 

X-semi axis of 
ellipsoid 

ft 

100 

B 

IDXV 

+ 

1 

91 

Y-semi axis of 
ellipsoid 

ft 

100 

C 

mxv 

+ 

2 

92 

Z-semi axis of 
ellipsoid 

ft 

1 

AP(l) 

IDXV 

+ 

3 

93 

S/S frame position 

ft 

0 

AP(2) 

IDXV 

+ 

h 

9U 

coordinates of antenna 

ft 

0 

AP(3) 

I 

IDXV 

+ 

5 

9$ 

gibmal point 
* 

ft 

- 50 

IR0T 

IDXV 

+ 

6 

96 

rotation sequence for 
following angles 

— 

321 

hii(i) 

IDXV 

+ 

7 

97 

ellipsoid X-axis 
rotation angle 

deg 

0 

hh(2) 

IDXV 

+ 

8 

98 

ellipsoid Y-axis 
rotation angle 

deg 

0 

HII(3) 

IDXV 

+ 

9 

99 

ellipsoid Z-axis 

deg 

0 


rotation angle 


A -21 



Input List For Subroutine 3ENTSS2 


Variable Relative Absolute Nominal 

Name Location Location Description Units Value 


A 

IDXV 



100 

X-semi axis of ellipsoid 

ft 

100 

B 

IDXV +1 

101 

Y-semi axis of ellipsoid 

ft 

100 

C 

IDXV 

+ 

A 

C. 

102 

Z-semi axis of ellipsoid 

ft 

1 

AP(l) 

IDXV 

+ 

3 

103 

S/S frame position coor- 

ft 

0 

AP(2) 

IDXV 

+ 

k 

10U 

dinates of antenna gim- 

ft 

0 

AP(3) 

IDXV 

+ 

5 

105 

bal point 

ft 

- 5o 

IR0T 

IDXV 

+ 

6 

106 

rotation sequence for 
following angles 

— 

321 

PHI(1) 

IDXV 

+ 

7 

10? 

ellipsoid X-axis rota- 
tion angle 

deg . 

0 

PHI(2) 

IDXV 

+ 

8 

108 

ellipsoid Y-axis rota- 
tion angle 

deg 

0 

PHI(3) 

IDXV 

+ 

9 

109 

ellipsoid Z-axis rota- 

deg 

0 


tion angle 
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Input List For Subroutine- ANTDXN 

Variable Relative Absolute Nominal 


Name 

Location 

Location 

Description 

Units 

Value 

IAA(1, 1) 

IBXV 



600 


slug ft^ 

1000 

IAA(2, 2) 

IDXV 

+ 

1 

601 


slug ft ^ 

1000 

IAA(3, 3) 

IDXV 

+ 

2 

602 

> inertia diadic of antenna 

slug ft^ 
0 

2000 

IAA(1, 2) 

IDXV 

+ 

3 

603 j 

in antenna reference 

slug ft 

0 

IAA(1, 3) 

IDXV 

+ 

U 

60U 

iframe 

slug ft 2 

0 

IA.A(2, 3) 

IDXV 

+ 


605 1 

I 

slug ft 2 

0 

MBAR 

IDXV 

+ 

6 

60 6 

reduced mass 

slugs 

10 

R2A(1) 

IDXV 

+ 

7 

607 ) 

1 antenna center of mass to 

ft 

0 

R2A(2) 

IDXV 

+ 

8 

608 

>hinge pt. vector in 

ft . 

0 

R2A(3) 

IDXV 

+ 

9 

609 ! 

(antenna reference frame 

ft 

10 

EBETAA(l) 

IDXV 

+ 

10 

610 'j 

IBETA gimbal axis direc- 

— 

0 

EBETAA(2) 

IDXV 

+ 

11 

611 

tion vector in antenna 

■ — 

1 

EBETAA(3) 

IDXV 

+ 

12 

612 | 

( reference frame 

- 

0 

RlP(l) 

IDXV 

+ 

13 

613 j 

1 S/S center of mass to 

ft 

0 

R1P( 2 ) 

IDXV 

+ 

1U 

6ia 

►hinge pt. vector in 

ft 

0 

R1P(3) 

IDXV 

+ 

15 

615 ) 

(S/S reference frame 

ft 

-SO 

EALFAP 

IDXV 

+ 

16 

616 ) 

^ALFA gimbal axis direc- 

— 

1 

EALFAP 

IDXV 

+ 

17 

617 

tion vector in S/S 

— 

0 

EALFAP 

IDXV 

+ 

Q 

18 

618 ' 

(reference frame 

— 

0 

ALFA 

IDXV 

+ 

19 

619 

initial ALFA gimbal angle 

deg 

- 8U.86 

ALFAD 

IDXV 

+ 

20 

620 

initial ALFA gimbal rate 

deg/sec 

0 

BETA 

IDXV 

+ 

21 

621 

initial BETA gimbal angle 

deg 

0 

BET AD 

IDXV 

+ 

22 

622 

initial BETA gimbal rate 

deg/sec 

0 

FREQ 

IDXV 

+ 

23 

623 

bending mode frequency 

rad/sec 

.628 


For FREQ ^ 0 no S/S 
flexibility is accounted 
for and no further inputs 
are required 
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Input List For Subroutine ANTDYN Continued 


Variable 

Relative 

Absolute 


Nominal 

Name 

Location 

Location 

1 Description 

Units 

Value 




For Freq > 0 



MP 

IDXV + 21* 

621* 

S/S mass 

slugs 

10000 

mrh 

IDXV + 2 5 

625 \ 


ft 2 

13860 

chhp(i) 

hjxv + 26 

626 


ft 

0 

CPHIP(2) 

3DXV + 27 

627 


ft 

- 1 

CH1IP(3) 

IDXV + 28 

628 

1 

ft 

0 

CYP(l) 

IDXV + 2 9 

629 


ft 2 

77.88 

CYP(2) 

IDXV + 30 

630 

• 

ft 2 

0 

CYP(3) 

IDXV + 31 

631 

t 

ft 2 

0 

CEP(1, 1) 

IDXV + 32 

632 

' 

ft 2 

.1*0 

CEP(1, 2) 

IDXV + 33 

633 


ft 2 

0 

CEP(1, 3) 

IDXV + 31* 

631* 


ft 2 

0 

CEP(2, 2) 

IDXV + 35 

635 

l flexible S/S mass and 

ft 2 

.223 

CEP(2, 3) 

IDXV + 36 

636 

\bending deflection dis- 
tribution parameters 

ft 2 

0 

CEP(3j 3) 

IDXV + 37 

. 637 

ft^ 

. 2 ia 

CNBP(1, 1) 

IDXV + 38 

638 

1 

ft 2 

0 

CNBP(1, 2) 

IDXV + 39 

639 


. ft 2 

0 

CNBP(1, 3) 

IDXV + 1*0 

61*0 


ft 2 

0 

CNBP(2, 1) 

IDXV + ia 

6ia 


ft 2 

0 

CNBP(2, 2) 

IDXV + 1*2 

61*2 


ft 2 

0 

CNBP(2, 3) 

IDXV + 1*3 

61*3 


ft 2 

62.76 

CNBP( 3 > 1) 

IDXV + 1*1* 

61*1* 


ft 2 

0 

CNBP(3, 2) 

IDXV + h$ 

61*5 I 

: 

ft 2 

12.76 

CNBP(3, 3) 

3D XV + 1*6 

61*6 / 


ft 2 

0 

PSICP(l) 

IDXV + 1*7 

61*7 ' 

bending deformation slope 

— 

.0175 

FSICP(2) 

IDXV + 1*8 

61*8 

► at control position 

— 

0 

PSICP(3) 

IDXV + 1*9 

! 61*9 J 


- 

0 
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Input List For Subroutine ANTDYN Continued 


Variable 

Relative 

Absolute 



Nominal 

Name 

Location 

Location 

Description 

Units 

Value 

FSIGP(l) 

IDXV + 50 

650 ) 

| 


.0175 

PSIGP(2) 

IDXV + 51 

651 

Lbending deformation slope 
[at IMU position 

— 

0 

PSIGP(3) 

IDXV + 52 

652 ! 

> 

— 

0 

Q 

IDXV + 53 

653 

initial bending coordinate 
value 

— 

0 

QD 

IDXV + Sk 

6$h 

initial bending coordinate 

-1 

sec 

0 


rate 
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Input List For Subroutine CNTRL7 




Variable Relative Absolute 

Name Location Location 

Description 

Units 

Nominal 

Value 


ALFA channel parameters: 


KLA 

3NDEXV 

iia 

attitude gain, pointing 
mode 

deg -1 

vri 

O 

K2A 

INDEXV + 1 

11*2 

rate gain, pointing mode 

(deg/sec) 1 

5oo 

K3A 

INDEXV + 2 

11*3 

attitude gain, tracking 
mode 

deg* 1 

5o 

KUA 

INDEXV + 3 

11*1* 

rate gain, tracking mode 

(deg/sec) -1 

5oo 

k5a 

INDEXV + h 
# 

11*5 

integral gain, tracking 
mode 

sec*" 1 

CsJ 

O 

• 

alfab 

INDEXV + 5 

ll*6 

1 

bias 

deg 

0 




BETA channel parameters: 



JOB 

INDEXV + 6 

11*7 

attitude gain, pointing 
mode 

deg' 1 

5o 

K2B 

INDEXV + 7 

11*8 

rate gain, pointing mode 

(deg/sec) -1 

5oo 

KJB 

INDEXV + 8 

11*9 

attitude gain, tracking 
mode 

deg" 1 

5o 

KUB 

INDEXV + 9 

1$0 

rate gain, tracking mode 

(deg/sec) 1 

5oo 

K5B 

INDEXV + 10 

151 

integral gain, tracking 
mode 

-1 

sec 

.02 

BETAB 

INDEXV + 11 

152 

bias 

deg 

0 




resolver characteristics: 



HES0LP 

INDEXV + 12 

153 

position quantization 

deg 

0 

RES0LR 

INDEXV + 13 

151* 

rate quantization 

deg/sec 

0 

RE SOLI 

INDEXV ♦ ll* 

155 

output signal quantiza- 
tion 


0 




spiral scan parameters: 



TSURCH 

INDEXV + 1 ! $ 

156 

allowed search duration 

sec 

5 

R0 

INDEXV + 16 

157 

initial scan radius 

deg 

.05 

VT 

INDEXV + 17 

158 

tangential scan rate 

deg/sec 

.8 

DS 

INDEXV + 18 

159 

radial scan rate 

deg/rev 

1.0 
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Input List For Subroutine 


CNTRL7 Continued 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 




initial conditions: 



REGA 

INDEXV + 19 

160 

ALFA channel integrator 

— 

0 

REGB 

INDEXV + 20 

161 

BETA channel integrator 

— 

0 

QUANTM 

INDEXV + 21 

162 

# of bits in gimbal sen- 
sor output 

— 

0 

QUANT P 

INDEXV + 22 

163 

§ of bits in computed 
gimbal output 


0 

QUANTZ 

INDEXV + 23 
» 

16 U 

# of bits in error signal 
output 


0 

LSTJRCH 

INDEXV + 2h 

165 

flag: 0 to inhibit spiral 
scan 

— 

0 
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Input List For Subroutine SERV0, FGEN 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 

KGA 

INDEXV 

181 

ALFA channel torque gain 

ft -lbs 

1 

KGB 

INDEXV + 1 

182 

BETA channel torque gain 

ft -lbs 

1 




Channel 1 of FGEN 



NFTS 

INDEXV + 2 

183 

• 

# of points 

— 

3 

XP(l,l) 

INDEXV +3 

18U 

input value, pt. 1 

rad/sec 

0 

YP(l,l) 

INDEXV + h 

185 

output value, pt. 1 

ft -lbs 

.16 

XP(1,2) 

INDEXV + 5 

186 

input value, pt. 2 

rad/sec 

•05 

YP(l,2) 

INDEXV + 6 

187 

output value, pt. 2 

ft-lbs 1 

CO 

O 

• 

XP(1,3) 

INDEXV + 7 

188 

input value, pt. 3 

rad/sec 

10. 

YP(1,3) 

INDEXV + 8 

189 

output value, pt. 3 

ft -lbs 

.08 




Channel 2 of FGEN 



NFTS 

INDEXV + 9 

190 

# of points 

- 

3 

XP(2,1) 

INDEXV + 10 

191 

input value, pt. 1 

rad/sec 

0 

YP(2,1) 

INDEXV + 11 

192 

output value, pt. 1 

ft-lbs 

.16 

XP(2,2) 

INDEXV + 12 

193 

input value, pt. 2 

rad/sec 

.05 

YP(2,2) 

INDEXV + 13 

194 

output value, pt. 2 

ft-lbs 

.08 

XP(2,3) 

INDEXV + 1U 

195 

input value, pt. 3 

rad/sec 

10. 

YP(2,3) 

INDEXV + 1 5 

196 

output value, pt. 3 

ft-lbs 

.08 
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Input List For Surboutine 0BKNMT 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 

IR0T 

INDEXV 

32$ 

input rotation sequence 

— 

321 

PHI(l) 

INDEXV + 1 

326 

S/S X-axis rotation angle 

deg 

-10. 

HII(2) 

INDEXV + 2 

327 

S/S Y-axis rotation angle 

deg 

0 

PHI(3) 

INDEXV + 3 
1 

oo 

CM 

S/S Z-axis> rotation angle 

* • 

t 

deg 

0 
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Input List For Subroutine 0EEPH3 


Variable 

Name 

Relative 

Location 


Absolute 

Location 

Description 

Units 

Nominal 

Value 

K0FT 

3NDXV 



500 

option index (not used) 

— 

— 

XT(2) 

INDEXV 

+ 

1 

501 

time at input point 1 

sec 

0 

RV1(2) 

INDEXV 

+ 

2 

502 


m.m. 

0 

RV2(2) 

INDEXV 

+ 

3 

503 

earth frame position 
coordinates of point 1 

m.m. 

38UU.7 

RV3(2) 

INDEXV 

+ 

h 

5oli 

m.m. 

0 

Wl(2) 

INDEXV 

+ 

$ 

505 


ft/sec 

0 

W2(2) 

INDEXV 

i 

+ 

6 

506 

|earth frame velocity . 
icomponents at point 1 

ft/sec 

0 

W3(2) 

INDEXV 

+ 

7 

507 

! i 

ft/sec 

21*552.3 

XT(3) 

INDEXV 

+ 

11 

5n 

time at input point 2 

sec 

1200. 

RV1(3) 

INDEXV 

+ 

12 

512 


m.m. 

102.1*1 

RV2(3) 

INDEXV 

+ 

13 

5i3 

earth frame position 
coordinates of point 2 

m.m. 

1167-31 

RV3(3) 

INDEXV 

+ 

lh 

5U* 

m.m. 

3660.26 

Wl(3) 

INDEXV 

+ 

15 

5i5 


ft/sec 

-201*2.1 

W2(3) 

INDEXV 

+ 

16 

5i6 

'earth frame velocity 
|C components at point 2 

ft/sec 

-23277.3 

W3(3) 

INDEXV 

+ 

17 

517 

ft/sec 

7U7U.1 

XT(1*) 

INDEXV 

+ 

21 

521 

time at input point 3 

sec 

21*00. 

RVl(U) 

INDEXV 

+ 

22 

522 


m.m. 

-51*1*. 82 

RV2(U) 

INDEXV 

+ 

23 

523 1 

earth frame position 
coordinates of point 3 

m.m. 

-3081.21 

RV3(U) 

INDEXV 

+ 

2h 

52U I 

m.m. 

2233.22 

Wl(l*) 

INDEXV 

+ 

25 

525 


ft/sec 

-2l*81*.2 

W2(l*> 

INDEXV 

+ 

26 

526 

earth frame velocity 
components of point 3 

ft /sec 

-11*01*9 . 3 

W3(ii) 

INDEXV 

+ 

27 

527 


ft/sec 

-19972.7 
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Input List For Subroutine 0EEFH3 Continued 

Variable Relative Absolute Nominal 


Name 

Location 

Location 

Description 

Units 

Value 

XT(5) 

INDXV 

+ 31 

531 

time at input point 1* 

sec 

3600. 

RVl(5) 

INDXV 

+ 32 

532 

| 

m.m. 

-799.1*1* 

RV2(5) 

INDXV 

+ 33 

533 

[earth frame position 
[coordinates of point 1* 

m.m 

-2975.0 

RV3(5) 

INDXV 

+ 3U 

531* 

[ » 

m.m. 

-2299.55 

wi(5) 

INDXV 

+ 35 

535 


ft/sec 

3812.1; 

W2(5) 

INDXV 

+ 36 

536 

[earth frame velocity 
| components at point .1* . 

ft /sec 

11*187-3 

W3(5) 

INDXV 

* 

+. 37 

537 


ft/sec 

-19662.8 

XT(6) 

INDXV 

+ I4I 

5 1 * 1 

\ 

time at input point 5 

sec 

U800. 

RV1(6) 

INDXV 

+ h2 

51*2 


m.m. 

U28.52 

RV2(6) 

INDXV 

+ k3 

51*3 

earth frame position 
coordinates of point 5 

m.m. 

1173.86 

RV3(6) 

INDXV 

+ .1*1*“ 

51*1* 

m.m. 

-3631 *.1*5 

Wl(6) 

INDXV 

+ 1*5 

51*5 


ft/sec 

7965-5 

W2(6) 

INDXV 

+ 1x6 

51*6 

earth frame velocity 
components at point 5 

ft/sec 

21795.1* 

W3(6) 

INDXV 

+ 1*7 

51*7 

ft/sec 

7970.7 

xt(7) 

INDXV 

+ 51 

551 

time at input point 6 

sec 

6000. 

RV1(7) 

INDXV 

+ 52 

552 

1 

m.m. 

1628.66 

RV2(7) 

INDXV 

+ 53 

553 

[earth frame position 
jcoordinates of point 6 

m.m. 

31*81.80 

RV3(7) 

INDXV 

+ 51* 

551* 

m.m. 

82.13 

Wl(7) 

INDXV 

+ 55 

555 

| 

ft/sec 

-222.1; 

W2(7) 

INDXV 

+ 56 

556 

[earth frame velocity 
[components at point 6 

ft/sec 

-1*75.5 

W3(7) 

INDXV 

+ 57 

557 

I 

ft/sec 

21*51*6.7 
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Input List For Surboutine 0EEPH1/SINCH 


Variable Relative Absolute 

Name Location Location 


Description 


Units 


K0FT 3NDXV 560 

RL0M INDXV + 1 561 


option index (not used) — 

relay satellite longitude deg 






Nominal 

Value 


{ 15 ' 


A -32 



Input List For Subroutine 0BEFH2/SXNCH 


Variable Relative Absolute 

Name Location Location 


Descripti ->n 


Units 


K0PT INDXV 580 option index (not used) — 

RL0N INDXV + 1 581 relay satellite longitude deg 


Nominal 

Value 
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Input List For Subroutine 


CLBRTN 


Variable 

Name 

Relative 

Location 

Absolute 

Location 

Description 

Units 

Nominal 

Value 

KfifPT 

IDXV 


3li0 

option index (not used) 

— 

— 

izMega 

IDXV 

+ 1 

3hl 

frequency of sinusoidal 
corrections 

deg/sec 

0 

tref 

IDXV 

+ 2 

3U2 

reference time 

sec 

0 

AA 

IDXV 

+ 3 

3i*3 

ALFA gimbal angle bias 
corrections 

deg 

0 

BA 

IDXV 

+ h 

3hh 

Coeff. of sine correc- 
tion term (ALFA) 

deg 

0 

CA 

IDXV 

+ 5 

3U5 

Coeff. of cosine correc- 
tion term (ALFA) 

deg 

0 

AB 

IDXV 

+ 6 

3U6 

BETA gimbal angle bias 
correction 

deg 

0 

BB 

raxv 

+ 7 

3U7 

Coeff. of sine correc- 
tion term (beta) 

deg 

0 

CB 

IDXV 

+ 8 

3U8 

Coeff. of cosine correc- 
tion term (BETA) 

deg 

0 




nw uvo 
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Figure A-l Functional Block Diagram of Setup/Loop 
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APPENDIX B 

BRIEF DESCRIPTION OF THE DEADAM PACKAGE 

A real-time simulation analysis can be accomplished with LMSC DEADAM simula- 
tion. Figure B-l illustrates the basic structure of the program. The main 
program DEADAM acts as the overall executive program which controls the main 
integration package (ADAMS) and the primary I/O routines. The ADAMS sub- 
routine contains the heart of the integration package. It performs numerical 
integration of the system of first order ordinary differential equations with 
piecewise continuous derivatives. The differential equations are usually non- 
linear and time varying and can be represented in the following form. 

X = f (x, u(t)) 

The above differential equations are defined by the user through the con- 
struction of the user executive (SETUP). The integration is performed using 
a third-order Adams -Bashforth predictor-corrector formula. The step size is 
adjusted automatically to stay within prescribed error bounds on the difference 
between predictor and corrector. Hence the step size will vary to allow 
efficient use of computer time without sacrificing accuracy. A self -starting 
procedure is used to form the first few differences. 

The above features are standard for most Adams -Bashforth methods. In addition, 
a discontinuity prediction scheme has been incorporated into the DEADAM pack- 
age. This allows the discontinuity point to be approached without iteration 
which results again in efficient use of compute^ time. 
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Figure B-l DEADAM Program 





APHSNDIX B 

The user executive (SETUP) must be written for each simulation by the user . 

It is the organizing subroutine through which the differential equations are 
defined. The prescription for writing SETUP is contained in the document: 

"DEADAM" A Program Especially Suited for Nonlinear Simulation 
of Attitude Control Systems," ASCS/2 Li 8£/6211, 29 Aug. 1969* 

Since many aerospace control system simulations have parts in common, a user 
library has been developed for rapid and easy construction of complicated 
simulations. Partial lists with brief descriptions of the existing library 
subroutines are contained in Tables B-2 and B-3. The standarization of these 
routines allows the construction of sophisticated simulations of highly com- 
plex systems in very short periods of time. 

The ease in construction of a typical DEADAM simulation can best be illustrated 
via a simple example. Suppose a simulation was required for. a rigid vehicle 
stabilized by a pitch and yaw wheel during a boost thrust along its. x axis. 

Also suppose there exists a pitch and yaw torque on the rigid body. The pro- 
blem is to observe the motion of the body and the change in the orbital states. 

The user may write a SETUP routine taking advantage of the user library. 

Such a SETUP is listed on, the next page . Note that other outputs which may 
be desired have not been defined . 

There are many fixed input variables which are not shown explicitly in the 
context of the user executive (SETUP). These are defined indirectly frcm an 
input array through the initialization calls to 0RBDYN and KINDYN. Examples 
of such variables are inertias, initial rigid body orientation, initial rates, 
inital orbit states. Also initialized through the same input array are the 
fixed parameters needed by the integration package and those needed by the 
1/0 routines. Examples of these are minimum step size, print interval, initial 
time, maximum time and number of output channels. 
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SUBROUTINE SETUP (KP) 

C<ZMMON/lO/OUT (£0, It 8 ) 

C0MM0N/C0NST/DEG, RAD, PI, TPI, PI2, REQ, RPL, GM, QMGE, FTNM 

cOmmOn/setOut/ichanl( 96 )/sOd/inpadm/main/aid( 13 ) , V( 300 ) 
C0MM0N/KNDN/TBO(3>3), JUNK3(18), PHI(3), RATES(3), 0MEGA(3) 
C0MM0N/0RBXT/WREF(3), RMAG, RV(3), t0e(9) vmag, w(3), alt, LAT, LON 
REAL LAT, L0N, T0RQB(3), HWHEEL(3), F0RCE(3), THRUST(3)> TE0(3,3) 

GO TO ( 100 , 200, 300 , 1*00, $ 00 , 600 , 700 ,), KP 


100 lO 

= 

1 

VIndex of initial state . variable . 

INPADM 

= 

100 VIndex for ADAMS input data 

RETURN 
200 NVORB 


1 

VIndex for ORBDYN input data 

NVKIN 

= 

50 

VIndex for KINDYN input data 

IDXORB 


1 

V0RBDYN state index 

CALL ORBDYN (1, 

NVORB, IDXORB, FORCE, SMUG) V Initialization 

IDXKIN 

= 

8 

VKINDYN state index 


CALL KINDYN ( 1 , NVKIN, IDXKIN, TORQB, HWHEEL, WRED, RMAG) V Initialization 

NND = 1U VIndex of final state variable 

RETURN 

300 CALL ORBDYN (2,NV0RB, IDXORB, FORCE, SMUG) V WRITE OUT 

CALL KINDYN (2, NVKIN, IDXKIN, T0RQB, HWHEEL, WREF, RMAG) V WRITE 0UT 
ICHANL(l)= 6H ROLL VAdditional output channels heading 
ICHANL( 2 ) = 6H RATE Vcan be defined here 
RETURN 

1*00 CALL 0REDYN(3, NVORB, IDXORB, FORCE, SMUG) V DEFINE OUTPUT 

CALL KINDYN(3, NVKIN, IDXKIN, TORQB, HWHEEL, WREF, RMAG) V DEFINE OUTPUT 
THRUST(l)= 200. VThrust in lb. 

SMUG = 100. VMass in slugs 

CALL TRNPS(TOB, TBO) 

CALL DDV(TBO, THRUST, F0RCE) 

CALL. ORBDYN (1*, NVORB, IDXORB, FORCE, SMUG) V DEFINE DERIVATIVES 
T0RQB(2) = .5 VPitch torque in ft -lb 

T0RQB(3) = .6 VYaw torque in ft-lb 

HWHEEL(2)= $0. VPitch wheel angular momentum ft-lb-sec 

HWHEEL(3)= 5>0. VYaw wheel angular momentum ft-lb-sec 

CALL KINDYN (1*, NVKIN, IDXKIN, T0RQB, HWEEL, WREF, RMAG) 

RETURN 

$00 CALL 0RBDYN(5, NVORB, IDXORB, FORCE, SMUG) V SPECIAL OUTPUTS 
CALL KINDYN($, NVKIN, IDXKIN, TORQB, HWHEEL, WREF, RMAG) 

OUT(LOUT,l) = RATES(l)*DEG VAdditional outputs defined here 
RETURN 
600 RETURN 

700 CALL 0RBDYN(6, NVORB, IDXORB, FORCE, SMUG) V DEBUG 

CALL KINDYN(6, NVKIN, IDXKIN, T0RQB, HWHEEL, WREF, RMAG) V DEBUG 
RETURN 

Once the user executive (SETUP is written and compiled, the initiation of the 
program is carried out by 
XQT DEADAM 

Following by appropriate 
DATA CARDS 
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Figure B-2 Functional Block Diagram of Example Simulation 
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APPENDIX B 

TABLE B-2 

USER LIBRARY - AUXILIARY SUBROUTINES (AND FUNCTIONS ) 


Subroutine Name 


Brief Description 


ANGLE (X,Y) 
ANGMOD (X) 
ARCTAN (A, B) 


Defines angle between X and Y 
X module 2i r in radians 
Arc tangent function with U quadrant return (ATAN2) 


ANGLET (A, IAX, D) D = [A 3 (transformation for rotation A about axis IAX) 

AXEQB$(A,X,B) Solves AX = B for X (5 x l) 

B0XFRD(X1,X2,X3) Function = XI, X2, X3 

DCTOEA (KS,LS,MS,E,A,B,C,IDEG) Direction cosing matrix (E) to 

DCTOEP(D,EP,DMAG) Direction cosing matrix to Euler parameters 
DDD(D1,D2,D3) D3 = D1*D2 Matrix product (see MAB) 

DDV(D,V1,V2) V2 = D .VI Linear transformation (3x3) of a vector 

DIAD(R,S,RST) RST = R > < S (Diadic) 

DINV (D,Dl) DI -'(D) 

DOTPRD(X,Y) Dot product of vectors X and Y 

EATODC(KS,LS,MS,DC,AN,EN,CN,IDEG) DC = [CNj^tBNl^LAN]^ ((DEG = 0 rad.) 

EFTODC(EP,D) Euler parameters to direction cosines 

EULER ( IS, JS,KS,WBODY, WREF,A,R,DC,IDEG,IEULER) 

IDEG = 0, converts Euler angles (A) and rates (R) to body 
rates (VIBODY) and dir. cos (DC) 

IDEG = 1, converts (WBODY) and (DC) to (A) and (R), order 
of A is IS ,JS,KS 

EAYRP(PSI, PHI, THETA, T) Direction cosing matrix (T) to Euler angles 

ce ] 2 [<p ] 1 [t ) 3 

EPRATE(EPD,EPW) Rigid body rate (W) to Euler parameter rate (EPD) 

EPWINI (Q,W,U,PSI, PHI, THETA, UD,F5 ID, PHID,THETAD) Define inertial euler 

parameters (q) and body rate (W) 

EP2TB1 (TBI,Q) Euler parameters (Q) to direction cosine matrix (TBl) 
IENGTH (X) Real function determines magnitude of vector X (3xl) 


MAC (C,A,B) 

C = A*B 

MABC (D,A,B.C) 

D = A * 

MABT (C ,A ,B) 

C = A * 

MALAT (C,A,l) 

c = a t i 

MATDA (C,D,A) 

c = a t b 


3 T 

.T 


C = A^B where B (I,J ) • = A(l, J) * D (I) 
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TABLE B-2 (Continued) 

MATIA (C,A,l) C = A T IA 

MFILL (C,C11, C12, C13, C21, C22, C23, C31, C32, C33) Fills (3x3) matrix 
MINV (C,A) C = (A) (see DINV) 

MSUM2 (D,SA,A,SB,D) D = SA * A '+ SB * B sums two matrices (SA 3 SB = scalors) 
MSIM2 (D,SA,A,SB,B,SC,C) T sums three matrices (SA,SB, SC = scalors) 

MTR (C,A) C = A 1 

PAGE Skips to new page 

TFILL (A, NAXIS > ANGLE) A = [Angle] NAXIg (see ANGLET) 

TF12 (A, ANGLE1, ANGLE 2) F0RMS Roll/Pitch or 1/2 rotation 
TRNPS (D,DT) DT = (D) t (see MTR) 

VABCTX (V,A,B,C,X) V = (A B C) T X 

VABCX (V,A,B,C,X) V = (A B C).X 

VABTX (V,A,B,,X) V = (AB)^X 
VABX (V,A,B,X) V = (AB).X. 

VATX (V,A,X) V = A T X 

VATXSY ( VjA^XjSjY) V = A T X + S.Y (S = scaler 

VAX (VjA,X) V = AX (see DDV 

VAXSY (V,A,X,S,Y) V = AX + S . Y (S = scaler 

VCOPY (V,I,J,T) Copies row (I = row) or col. (I = col.) of T into V 

VCROSS (ZjS,XjY) Z = S . (X cross Y) cross product (S = scalor) (See VXV) 

VDX (V, D,X) V(l) = D(l) . X(I) I = 1,3 

VFILL (Z,Z1,Z2,Z3) Fills (s x l) vector 

VMAG (X) Forms magnitude of X (function) 

VSCALE (X,S) Rescales X to magnitude S 

VSUM2 (Z,S1,X1,S2,X2) Z = SI . XI + Se . X2 (SI, S2 = scalors) sums two vectors 
VSUM3 (S > S1,X1,S2,X2,S3,X3) (see VSUM2) 

VSUMli (S,S1 J X1,S2,X2,S3 J X3,SU,X1;) (see VSUM2) 

VUNIT (X) Rescales X to unit vector 

VXKYSZ( V ,X , Y,S , Z ) V = X cross Y + S.Z (S = scalor) 

VXDDV(V1,D,V2V3) V3 = VI cross D.V2 

VXV (Vi,W2,V3) V3 = VI cross V2 (cross product) (see VCROSS) 
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TABLE B-3 - S FECIAL OUTPUT ROUTINES 

MLOOK1 (T1,A) 

ML00K2 (T1,A1,T2,A2) 

MLOOK3 (T1,A1,T2,A2,T3,A3) 

MLO.QK1* (Tl,Al,T2,A2,T3,A3,Tli,AU) 


Prints out (3x3) matrices 
(A^jAgjA^A^) with associated 

titles (T1,T2,T3,TU), titles 
. are 30 spaces long, numbers are 
in FIO .3 format 


TLOOK1 (T1,A) 


TL00K2 (T1,A1,T2,A2) 


TLOOK3 (T1,A1,T2,A2,T3,A3) 
TLOOKit (T1,A1,T2,A2,T3,A3,TU,AU) 

Same as MLOOK series except 
F10.6 format 

VLOOK1 (T1,V1,I1) 

VLOOK2 (T1,V1,I1,T2,V2,I2) 

VLOOK3 (T1,V1,I1,T2,V2,I2,T3,V3,I3) 

VLOOKli (Tl,Vl,H,T2,V2,I2,T3,V3,I3,Tli,VU,lU) 

Prints out (3xl) vectors 
with associated title (2k space 
long) format is FI 3.6 times 1.0 
to the 11 power etc. (scaling) 

VLOOKN (T,V,N) 

Prints out (N x I) vector with 
associated title (2h spaces 
long) 

BITS (N, DELTA, X) 

Function converting X to n bit 
word with resolution DELTA 

SCALE (N, DELTS, X) 

Same as BITS but output is 
(0 - 100 ) percent of full 
scale ! 

TRIGGR (D1,D2,D3,DU,D5,D6,D7,D8) 

Function which converts 8 real 
variable discretes ( 0 ,l) to 3 
octal bit word. 
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APPENDIX C 


TRADE-OFF STUDY SIMULATION RESULTS 

During the course of developing and debugging the simulation program des- 
cribed in this report numerous runs were made using typical values for the 
input parameters of the simulated system. The parameters which determine 
the response of the antenna pointing control system were adjusted in order 
to obtain stable, non-oscillatory, and rapidly converging antenna pointing 
control. 

The plots and printed computer output included in this appendix show the 
behavior of a typical space station computer steered high gain antenna using 
"optimized" control system gains. The values for all of the input parameters 
used in the illustrated simulation run are included in the input lists of 
Appendix A. 

The initial conditions were chosen to give a pointing error of about 10° in 
order to demonstrate the behavior of the controlled antenna during acquisition. 
When the gimbal angles get within 1.0° of the on board computed values the 
control system switches to its auto-track mode, using the monopulse error 
signals (at about 18 seconds in this case). Using an antenna with about a 1.0° 
half power beam width and monopulse horn separations of 1.0° (total) the 
antenna control system was able to track the relay satellite within about 0.3°. 

The bending input data used in this run was computed from the expressions 
derived in Section 3«U«2.3? for a space station main body of cylindrical mass 
distribution ($0 ft radius x 100 ft length) weighing 322,000 lbs (10,000 slugs) 
and with a boom antenna mount which deflects 1.0° angularly for each 1.0 ft 
translationally. The bending frequency was set at 1A0 Hz (.628 rad/sec) which 
is considered fairly low (corresponding to a not very stiff space station and 
boom). 
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The described bending mode is observed (Fig. C.ii) to respond as expected to 
glmbal torques and space station orbital rate and to have insignificant 
(10 ^ deg) effect on pointing. 
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Figure C.2 Gimbal Rate (ALFAD-), Gimbal Torque (TALFA), and 
Monopulse Error Signal (VEL) v.s. Time 
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Figure C.l* Bending Mode Coordinate v.s . Derivative 
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The following six pages are repro- 
duced from computer printout and show 
details of the data plotted in Figures 
C-l through C-U- 
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